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Abstract 

We consider random d-regular graphs on N vertices, with degree d at least (loglV)^. We 
prove that the Green’s function of the adjacency matrix and the Stieltjes transform of its 
empirical spectral measure are well approximated by Wigner’s semicircle law, down to the 
optimal scale given by the typical eigenvalue spacing (up to a logarithmic correction). Aside 
from well-known consequences for the local eigenvalue distribution, this result implies the 
complete (isotropic) delocalization of all eigenvectors and a probabilistic version of quantum 
unique ergodicity. 


1. Introduction and results 


1.1. Introduction. Let A be the adjacency matrix of a random d-regular graph on N vertices. 
For fixed d ^ 3, it is well known that as A —>• oo the empirical spectral measure of A converges 
weakly to the Kesten-McKay law [301135], with density 

(1.1) 


Thus, the rescaled adjacency matrix (d — 1) has asymptotic spectral density 


Qd{x) 



1 

d- 1 


^ ^[4 - 3;^]+ 
d ) 27r 


( 1 . 2 ) 


Clearly, gd{x) —)• g{x) as d —)• 00 , where g{x) ■= ^\/[4 — is the density of Wigner’s semi¬ 
circle law. The semicircle law is the asymptotic eigenvalue distribution of a random Hermitian 
matrix with independent (upper-triangular) entries (correctly normalized and subject to mild 
tail assumptions). From (II.2p it is natural to expect that, for sequences of random d-regular 
graphs such that d —)• 00 as A —00 simultaneously, the spectral density of (d — A con¬ 

verges to the semicircle law. This was only proved recently [H] (in m it was also shown with 
the restriction that d is only permitted to grow logarithmically in A). 

In the study of universality of random matrix statistics, local versions of the semicircle law 
and its generalizations have played a crucial role; see for instance the survey [22]. The local 
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semicircle law is a far-reaching generalization of the weak convergence to the semicircle law 
mentioned above. First, the local law admits test functions whose support decreases with N 
so that far fewer than N eigenvalues are counted, ideally only slightly more than order 1. (In 
contrast, weak convergence of probability measures applies only to macroscopic test functions 
counting an order N eigenvalues). Second, the local law controls individual matrix entries of 
the Green’s function. Both of these improvements have proved of fundamental importance for 
applications. In particular, the local law established in this paper is a crucial input in [3], 
where, with J. Huang, we prove that the local eigenvalue statistics of A coincide with those 
of the Gaussian Orthogonal Ensemble; see also Section 11.41 below. For Wigner matrices, i.e. 
Hermitian random matrices with independent identically distributed upper-triangular entries, 
the semicircle law is known to hold down to the optimal spectral scale 1 /N, corresponding to the 
typical eigenvalue spacing, up to a logarithmic correction. In [2l l 1711^144] , it was shown that the 
semicircle law (for d —)• oo) or the Kesten-McKay law (for fixed d) holds for random d-regular 
graphs on spectral scales that are slightly smaller than the macroscopic scale 1 (typically by a 
logarithmic factor; see Section [1.41 below for more details). 

In this paper we show that d-regular graphs with degree d at least (log N)'^ obey the semicircle 
law down to spectral scales {logN)‘^/N. This scale is optimal up to the power of the logarithm. 

From the perspective of random matrix theory, the adjacency matrix of a random d-regular 
graph is a symmetric random matrix with nonnegative integer entries constrained so that all 
row and column sums are equal to d. These constraints impose nontrivial dependencies among 
the entries. For example, if the sum of the first k entries of a given row is d, the remaining 
entries of that row must be zero. Previous approaches to bypass this difficulty include local 
approximation of the random regular graph by a regular tree (for small degrees) and coupling 
to an Erdos-Renyi graph (for large degrees). These approaches have been shown to be effective 
for the study of several combinatorial properties, as well as global spectral properties of random 
regular graphs. However, they encounter serious difficulties when applied to the eigenvalue 
distribution on small scales (see Section oi below for more details). Our strategy instead 
relies on a multiscale iteration of a self-consistent equation, in part inspired by the approach 
for random matrices with independent entries initiated in [2T] and significantly improved in a 
sequence of subsequent papers (again see Section [L4l for details). In previous works on local 
laws for random matrices, independence of the matrix entries plays a crucial role in deriving 
the self-consistent equation (see e.g. m for a detailed account). While the independence of the 
matrix entries can presumably be replaced by weak or short-range dependence, the dependence 
structure of the entries of random regular graphs is global. Thus, instead of independence, our 
approach uses the well known invariance of the random regular graph under a dynamics of local 
switchings, via a local resampling of vertex neighbourhoods. We believe that our strategy of 
local resampling, using invariance under a local dynamics combined with a multiscale iteration, 
is generally applicable to the study of the local eigenvalue distribution of random matrix models 
with constraints. 

Notation. We use a = 0(b) to mean that there exists an absolute constant C > 0 such that 
|a| ^ Cb, and a S> 6 to mean that a ^ Cb for some sufficiently large absolute constant C > 0. 
Moreover, we abbreviate |a, 6] := [a,b] n Z. We use the standard notations a f\b ■■= min{a, 6} 
and aWb := max{a, b}. Every quantity that is not explicitly a constant may depend on N, which 
we almost always omit from our notation. Throughout the paper, we tacitly assume N ^ 1. 

1.2. Random regular graphs. We establish the local law for the following three standard models 
of random d-regular graphs. 
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Uniform model. Let N and d be positive integers such that Nd is even. The uniform model 
is the uniform probability measure on the set of all simple d-regular graphs on |l,iV]. (Here, 
simple means that the graph has no loops or multiple edges.) Equivalently, its adjacency matrix 
A is uniformly distributed over the symmetric matrices with entries in {0,1} such that all rows 
have sum d and the diagonal entries are zero. 

Permutation model. Let iV be a positive integer and d an even positive integer. Let ui,..., < 7^/2 
be independent uniformly distributed permutations on Sat, the symmetric group of order N. The 
permutation model is the random graph on N vertices obtained by adding an edge {i, o'^{i)} for 
each i G |1, A^l and // G |l,d/2]. Its adjacency matrix A is given by 

d/2 d 

■= Hj = 0-/.(*)) , (1.3) 

{ 1=1 { 1=1 


with the convention that ad-^ ~ for (i/2 + 1 ^ /r ^ d in the second equality. All vertices 
have even degree, and in general the graph may have loops as well as multiple edges. Each loop 
contributes two to the degree of its incident vertex. 

Matching model. Let N be an even positive integer and d a positive integer. Let ai,...,ad 
be independent uniformly distributed perfect matchings on |1,AA]. A perfect matching can be 
identified with a permutation of S'at whose cycles all have length two. As in the permutation 
model, a graph on |1, A^] is obtained by adding an edge {i, (T^{i)} for all i G |1, A"] and fi G |1, dj. 
Thus, the corresponding adjacency matrix is again 

d 

■= =^^*(0)- (1-4) 

/i=l 


Graphs of this model can have multiple edges but no loops. Their degree d is arbitrary, but 
their number of vertices must be even. 

The models introduced above include simple graphs (uniform model), graphs with loops and 
multiple edges (permutation model), and graphs with multiple edges but no loops (matching 
model). Throughout this paper, all statements apply to any of the above three models, unless 
explicitly stated otherwise. As discussed in Section 11.41 below, our approach is quite general, 
and applies to other models of random regular graphs as well. Eor brevity, however, we give the 
details for the three representative models introduced above. 

We shall give error bounds depending on the parameter 


A2 

D := d A —pr 

(uniform model), 

(1.5) 

A2 

D := dA — 
d 

(permutation and matching models) . 

(1.6) 


In particular, for the uniform model, D = d if d ^ \/A, and for the permutation and matching 
models, D = d if d ^ N. Throughout the paper, we make the tacit assumption D ^ 1, which 
leads to the conditions d ^ for the uniform model and d ^ for the permutation and 
matching models. 
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1.3. Main result. To state our main result, we first observe that the adjacency matrix A of any 
d-regular graph on N vertices has the eigenvector e := , 1)* with eigenvalue d, and 

that (by the Perron-Frobenius theorem) all other eigenvalues are at most d in absolute value. 
The largest eigenvalue d of the eigenvector e is typically far from the other eigenvalues, and it 
is therefore convenient to set it to be zero. In addition, we rescale the adjacency matrix so that 
its eigenvalues are typically of order one. Hence, instead of A we consider 

H ■.= {d-l)-^/‘^{A-dee*) . (1.7) 


Clearly, A and H have the same eigenvectors, and the spectra of {d — 1)“^/^T and H coincide 
on the subspace orthogonal to e. 

Our main result is stated in terms of the Green’s function (or the resolvent) of i7, defined 
by 

G{z) := (1.8) 

for z £ C+. Here C+ := {E + ir/ : i? G R, r; > 0} denotes the upper half-plane. We always use 
the notation z = E + \r] for the real and imaginary parts of z G C+, and regard E = E(z) and 
r] = r]{z) as functions of z. 

For z G C+, let 

„(.) := /TTd, = + 

J X z Z 

be the Stieltjes transform of the semicircle law. Here the square root is chosen so that m(z) G C+ 
for ^ G C+, or, equivalently, to have a branch cut [—2, 2] and to satisfy Vz'^ — A ~ 2 ; as \z\ —>■ 00 . 
We shall control the errors using the parameter 




1 , 1 


( 1 . 10 ) 


and the function 


E,ir) := 


1 + 




A v^. 


( 1 . 11 ) 


where r G [0,1]. Away from the two edges 2 : = ±2 of the support of the semicircle law, i.e. for 
[z ± 2| ^ e for some e > 0, the function F is linearly bounded: Fz{r) = Oe(r). Near the edges, 
Fz{f) ^ provides a weaker bound. 

We now state our main result. 


Theorem 1.1 (Local semicircle law). Let G{z) be the Green’s function (II.8p of any of the 
models of random d-regular graphs introduced in Section UE^ Let ^ log^ S> (log A^)^ and 
Then, with probability at least 1 — 

m.&^\Gii{z) - m{z)\ = 0{Fz{f^{z))), max|Gij( 2 ;)| = 0 (^$( 2 ;)), (1-12) 

I 


simultaneously for all z G C+ such that rj S> 


The condition H ^ in the statement of Theorem 11.11 implies the following restrictions on 
the degree of the graphs: 


f d 
< d < 



(uniform model), 


(permutation and matching models). 


(1.13) 

(1.14) 
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Thus, for the smallest possible degree d and the smallest spectral scale rj for which Theo¬ 
rem 11.11 applies, the parameter ^ needs to be chosen as small as permitted, which is slightly 
smaller than (logA^)^. In particular, the local semicircle law holds for all rj ^ (logA^)^/A^ and 
all d ^ (logA^)^ satisfying d ^ A^^/^(log for the uniform model and d ^ N‘^{\ogN)~^ for 

the permutation and matching models. 

The estimates (|1.12p have a number of well-known consequences for the eigenvalues and 
eigenvectors of H, and hence also for those of A. Some of these are discussed below. In fact, by 
the exchangeability of random regular graphs. Theorem 11.11 actually implies an isotropic version 
of the local semicircle law, as well as corresponding isotropic versions of its consequences for the 
eigenvectors, such as isotropic delocalization and a probabilistic version of local quantum unique 
ergodicity. We discuss the isotropic modifications in Section [HI and restrict ourselves here to the 
standard basis of 

For instance. Theorem 11.11 implies that all eigenvectors are completely delocalized. 

Corollary 1.2 (Eigenvector delocalization). Under the assumptions of Theorem \l.ll 
with probability at least 1 — all -normalized eigenvectors of A or H have i°°-norm of 

size 0(^/\/]V). 

Proof. Since A and H have the same eigenvectors, it suffices to consider H. Let = {va,i)^=i, 
a = 1,... ,N denote an orthonormal basis of eigenvectors with Hv^ = AqV„. Let ^ be as in 
Theorem ll.il and set rj := C^‘^/N for some large enough constant C. Note that 

2 2 

T. < E (A, - A?)2 + = -|ImG«(A„ + i„), 

By Theorem 11.11 there exists an event of probability at least 1 — e~^ ^ such that for all i and 

a the right-hand side above is bounded by 

r]lmGii{Xa-\-iri) ^ ? 7 |m(A„ -\-ir])\+ 0{r]y^ ?$(Aq, -F ir])) ^ 2??, 


where we used the bound 

\m{z)\ ^ 1, (1.15) 

which follows easily from ()1.9[) . Thus v‘^^^2rj = /N) as claimed, concluding the proof. □ 

Next, Theorem 11.11 yields a semicircle law on small scales for the empirical spectral measure 
of H. The Stieltjes transform of the empirical spectral measure of H is defined by 


s{z) := 


1 


^ 1 


N ^ Xa- Z 
a=l 


1 

N 


N 




2 = 1 


where Xi,... ,Xi\f are the eigenvalues of H. Theorem 11.11 implies that 

s(z) = m{z) + 0{F^{^^{z))) 


(1.16) 


(1.17) 


with probability at least 1 — e Following a standard application of the Helffer-Sjostrand 

functional calculus along the lines of m Section 8.1], the following result may be deduced from 
(fLlTl) . 
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Corollary 1.3 (Semicircle law on small scales). Let 

f 1 ^ 

Qil) ■= / Q{x)dx, v{I) := 


denote the semicircle and empirical spectral measures, respectively, applied to an interval I. Fix 
a constant K > 0. Then, under the assumptions of Theorem \l.ll for any interval I C [—K,K] 
we have 


u{I) - g{I) 


O 




\I\ 

+ \i\ 





(1.18) 


with probability at least 1 —e where |/| denotes the length of I and k{I) ■= dist(/, {—2, 2}) 

the distance from I to the spectral edges ±2. 


Corollary 11.31 savs in particular that, in the bulk spectrum, the empirical spectral density of 
H is well approximated by the semicircle law down to spectral scales /N. Indeed, fix e > 0 
and suppose that I C [—2 + e, 2 —e], so that k(/) ^ e. Then the right-hand side of (11.181) is much 
smaller than q{I) provided that |/| i:|> We deduce that the distribution of the eigenvalues 

of H is very regular all the way down to the microscopic scale. Moreover, clumps of eigenvalues 
containing more than (logA^)^ eigenvalues are ruled out with high probability: any interval of 
length at most (\ogN)^/N contains with high probability at most 0((logiV)^) eigenvalues. 

Remark 1.4. The estimate (11.181) deteriorates near the edges, when k{I) is small. Here we do 
not aim for an optimal edge behaviour, and (11.181) can in fact be improved near the edges by a 
more refined application of (jl.l7p . For example, from (I1.17|) we also obtain the estimate 


u{i) - e{i) 


o 


Vl\i\ 


1 1 \ e' 

T»V4 + (Ar|i|)i/4j + N 


(1.19) 


with probability at least 1 — which is stronger than (11.181) when |/| and k{I) are small. 

Moreover, as explained in Remark 1 1.6 1 b elow. (I1.18P itself, and hence estimates of the form (11.191) . 
can be improved near the edges. We do not pursue these improvements here. 

Remark 1.5. Theorem 11.11 has a simple extension in which the condition rj S> f'^/N is dropped. 
Indeed, using Lemma l2.II below, it is easy to conclude that, under the assumptions of Theorem 
11.11 for any G C+ with rj = 0{f,‘^/N) we have the estimate \Gij{z) — Sijm(z)\ = 0(;^) with 
probability at least 1 — 


Remark 1.6. Up to the logarithmic correction we expect that the estimates (11.121) cannot 
be improved in the bulk of the support of the semicircle law, i.e. for |Fi| ^ 2 — e. On the other 
hand, (I1.12p is not optimal for |Fi| ^ 2 — e. For example, a simple extension of our proof allows 
one to show that the term ^<h(z) on the right-hand sides of ()1.12l) can be replaced by the smaller 
bound 


llmm{z) f 


Nr] 


+ 


Vd 


+ 


/^y/3 

\Nr]) 


( 1 . 20 ) 


In order to focus on the main ideas of this paper, we give the proof of the simpler estimate (11.121) . 
In Appendix!^ we sketch the required changes to obtain the improved error bound (jl.20p . The 
bound ()1.12p is sufficient for most applications, including Corollaries 1 1. 21[L3l Finally, we remark 
that all of our error bounds are designed with the regime of bounded z in mind; as z ^ oo, 
much better bounds can be easily obtained. We do not pursue this direction here. 
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1.4. Related results. We conclude this section with a discussion of some related results. The 
convergence of the empirical spectral measure of a random d-regular graph has been previously 
established on spectral scales slightly smaller than the macroscopic scale 1. More precisely, 
in m Theorem 1.6], the semicircle law is established down to the spectral scale d for 
d —)• oo. In [m Theorem 2 and Remark 1], the semicircle law is established down to the spectral 
scale (logA^)“^ for d = (logA^)"^ with 7 > 1, and the spectral scale 1/d for d = (logA^)"^ with 
7 < 1. In [21 Theorem 5.1], it is shown that for fixed d the Kesten-McKay law holds down to 
the spectral scale (logA^)”'^ for some c > 0. Finally, in [261 Theorem 2.1], it is shown that for 
fixed d the Kesten-McKay law holds down to the spectral scale (log 

The results of jU] were proved by coupling to an Erdds-Renyi graph. The probability that 
an Erdds-Renyi graph in which each edge is chosen independently with probability p is d-regular, 
with d = pN, is at least exp(—cA^logd). Hence, any statement that holds for the Erdds-Renyi 
graph with probability greater than 1 — exp(—cAIlogd) also holds for the random d-regular 
graph. While global spectral properties can be established with such high probabilities, super¬ 
exponential error probabilities are not expected to hold for local spectral properties. 

In a related direction, contiguity results imply that almost sure asymptotic properties of 
various models of random regular graphs can be related to each other (see e.g. [45] for details). 
Such results are difficult to extend to the case where d grows with N, for example because the 
probability that a graph of the permutation model is simple tends to zero roughly like exp(—cd^). 
This probability is smaller than the error probabilities that we establish in this paper. Our proof 
does not rely on a comparison between different models, but works directly with each model. It 
is rather general, and may in particular be adapted to other models of random regular graphs. 
Eor instance, by an argument similar to (but somewhat simpler than) the one given in Section [ 6 l 
we may prove Theorem II.II for the configuration model of random regular graphs. Moreover, by 
a straightforward extension of our method, our results remain valid for arbitrary superpositions 
of the models from Section 01 Eor example, we can consider a regular graph defined as the 
union of several independent uniform regular graphs of lower degree. (In fact, the matching 
model is the union of d independent copies of a uniform 1 -regular graph). 

The results of [2l ll7l[26] were obtained by local approximation by a tree. It is well known 
that, locally around almost all vertices, a random d-regular graph is well approximated by the 
d-regular tree, at least for fixed d ^ 3. The Kesten-McKay law is the spectral measure of the 
infinite d-regular tree, and many previous results on the spectral properties of d-regular graphs 
use some form of local approximation by the d-regular tree. In particular, it is known that the 
spectral measure of any sequence of graphs converging locally to the d-regular tree converges 
to the Kesten-McKay law; see for instance [9]. Moreover, in m, under an assumption on 
the number of small cycles (corresponding approximately to a locally tree-like structure and 
satisfied with high probability by random regular graphs), it is proved that eigenvectors cannot 
be localized in the following sense: if for some ^^-normalized eigenvector v = (uj)(^^ a set 
B C [l,iV] satisfies ^ e > 0, then \B\ ^ with high probability for some small 

5 (X e^. In comparision, for a random d-regular graph with d ^ (logA^)"^, Corollary 11.21 implies 
that if a set B has ^^-mass e > 0 then |H| ^ £:A^(log N)~'^ with high probability, which is optimal 
up to the power of the logarithmic correction. Eurthermore, in [2], for d-regular expander graphs 
with local tree structure, for fixed d ^ 3, a graph version of the quantum ergodicity theorem is 
proved; it is shown that averages over eigenvectors whose eigenvalues lie in an interval containing 
at least A^(log A^)“^ eigenvalues converge to the uniform distribution, along with a version of the 
Kesten-McKay law at spectral scales slightly smaller (by a related logarithmic factor) than the 
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macroscopic scale 1. For random regular graphs, also using the local tree approximation, similar 
estimates for eigenvalues on scales of order roughly (logA^)“^ were also established in |17[l26j . 
In all of these works, the logarithmic factor arises as the radius of the largest neighbourhood 
where the tree approximation holds, which is of order log^A^. 

Our proof does not use the tree approximation. Instead, we use that a local resampling using 
appropriately chosen switchings leaves the random regular graphs from Section 11.21 invariant. 
Switchings of random regular graphs were introduced to prove enumeration results in |36j : 
see also [45] for a survey of subsequent developments. Switchings are also commonly used 
for simulating random regular graphs using Monte Carlo methods; see e.g. m and references 
therein. Recently, switchings were employed to bound the singularity probability of directed 
random regular graphs [14| . 

For d-regular graphs, the value of second largest eigenvalue A 2 is of particular interest. At 
least for fixed d ^ 3, it was conjectured that for almost all random d-regular graphs we have 
A 2 = 2y/d — 1 + 0 ( 1 ) with high probability [T]. For fixed d, this conjecture was proved in [24] . 
following several larger bounds (for which references are given in [24] ). Very recently, the results 
of |24j were generalized and their proofs simplified in [71[39]. For the permutation model with 
d —)■ 00 as V —)■ 00 , the best known bound is A 2 = 0{y/d) [25] (see m Theorem 2.4] for a more 
detailed proof). 

Finally, it is believed that the eigenvalues of random d-regular graphs obey random matrix 
statistics as soon as d ^ 3. There is numerical evidence that the local spectral statistics in the 
bulk of the spectrum are governed by those of the Gaussian Orthogonal Ensemble (GOE) [28p38j . 
and further that the distribution of the appropriately rescaled second largest eigenvalue A 2 
converges to the Tracy-Widom distribution of the GOE [37] . 

In [3], with J. Huang, we prove that GOE eigenvalue statistics hold in the bulk for the 
uniform random d-regular graph with degree d £ [V", for arbitrary a > 0. Here, the 

lower bound on the degree is of purely technical nature, and we believe that the results of [3] 
can be established with the same method under the weaker assumption d ^ (logV)*^^^). The 
local law proved in this paper, in addition to the results of [271[34], is an essential input for the 
proof in [3] . 

Eor Erdds-Renyi graphs, in which each edge is chosen independently with probability p, the 
local semicircle was established under the condition piV ^ (log in [20]. Moreover, random 

matrix statistics for both the bulk eigenvalues and the second largest eigenvalue were established 
in [18] under the condition pN ^ iV^/s-i-a £qj. q- > 0 . Eor random matrix statistics of 

the bulk eigenvalues, the lower bound on pN was recently extended to pN ^ N°‘ for any a > 0 
in [27], and GOE statistics for the eigenvalue gaps was also established. Previous results on the 
spectral statistics of Erdds-Renyi graphs are discussed in [T81I201I44] . 


2. Preliminaries and the self-improving estimate 


In this section we introduce some basic tools and definitions on which our proof relies, and state 
a self-improving estimate. Proposition 12.21 from which Theorem 1 1.1 1 will easily follow. The rest 
of this paper will be devoted to the proof of Proposition 12.21 

Erom now on we frequently omit the spectral parameter z from our notation, and write 
G = G{z) and so on. The spectral representation of G implies the trivial bound 


Gij I ^ 


1 

V ' 


( 2 . 1 ) 





We shall also use the resolvent identity, for invertible matrices A,B 


= A-^{B - A)B-^. 


( 2 . 2 ) 


In particular, applying (|2.2I) to G — G*, we obtain the Ward identity 


EiG 


12 
ik \ 


Im Gii 


V 


(2.3) 


Assuming rj (|2.3p shows that the squared ^^-norm ^ '^k=i is smaller by the factor 

^ <C 1 than the diagonal element | Gu \. This identity was first used systematically in the proof 
of the local semicircle law for random matrices in |23] . 

The core of the proof is an induction on the spectral scale, where information about G is 
passed on from the scale r] to the scale r]/2. (See Remark 12.31 below for a comparison of this 
induction with the bootstrapping/continuity arguments used in the proofs of local laws in models 
with independent entries.) The next lemma is a simple deterministic result that allows us to 
propagate bounds on the Green’s function on a certain scale to weaker bounds on a smaller 
scale. This result will play a crucial role in the induction step. In order to state it, we introduce 
the random error parameters 


T = T{z) := max|Gij( 2 )| VI, T* = r*(z) := supT{E + irj'). (2.4) 

Lemma 2.1. For any M > 1 and z € C+ we have V{E + irj/M) ^ MT{E + ir]). 

Proof. Fix Fi G M and write r(r/) = r(£' + ir/). For sufficiently small h, since |x V 1 — y V 1| ^ 
|x — y| for X,y > 0, using the resolvent identity, the Cauchy-Schwarz inequality, and (12.3p . we 
get 


|F(t/+/ i) - F(7/)| ^ max\Gij{E+ i{ri + h)) - Gij{E+ ir])\ 




|/i|max^|Gjfc(F; + i(r/ + /i))Gfcj(F; + ir/)| ^ 


T{ri + h)r{rj) 


+ h)rj 


Thus, F is locally Lipschitz continuous, and its almost everywhere defined derivative satisfies 


aF 


drj 


F 

^ - 


This implies ^(yF(r/)) ^ 0 and therefore F(r//M) ^ MF(t/) as claimed. 


□ 


The main ingredient of the proof of Theorem 1 1.1 1 is the following result, whose proof consti¬ 
tutes the remainder of the paper. To state it, we introduce the set 


D = D(^) := ^ E + iri ^ rj ^ N , —N ^ E ^ N 


(2.5) 


where the implicit absolute constant in <C is chosen large enough in the proof of the following 
result. 
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Proposition 2.2. Suppose that ^ > 0, ( > 0, and that D ^ If for a fixed z gT) we have 

T*{z) = 0(1) 

with probability at least 1 — then for the same z we have 

max.\Gii{z) - m{z)\ = 0{F^{f,^{z))), m.ax\Gij{z)\ = 0(^$(z)) , (2.6) 

I i^j 

with probability at least 1 — e~(i^°s^)^C+ 0 {\ogN) ^ 

Given Proposition 12.21 Theorem 11.11 is a simple consequence. 

Proof of Theorem 11.11 Let ^log^ S> (log A^)^. We first note that it suffices to prove (I1.12p 
for z € D. Indeed, since d <C N'^ by assumption, the spectrum of H is contained in the interval 
[—d/{2^/df^^), dj (2V d — 1)] C [—and hence (jl.l2|) holds deterministically for \E\ ^ N 
by the spectral representation of G. Similarly, the proof of (I1.12p is trivial for ^ N. Since G 
is Lipschitz continuous in z with Lipschitz constant bounded by l/rf ^ A^^, it moreover suffices 
to prove ()1.12p for 2 e D n {N ^7?). By a union bound, it suffices to prove ()1.12l) for each 
E G [-N, N] n Ai-^Z. 

Fix therefore E G [-A^, A^] n (Ai-^Z). Let K := max{/c G N : A^/2^ ^ Cff^/N], where C > 0 
is the implicit absolute constant from the assumption p /N in the statement of the theorem. 
Clearly, K ^ 41og A^. For k G |0, AT], set pk •= N/‘2^ and z^ ■= E + ipf^. By induction on k, we 
shall prove that 

F*(zfc) ^ 2 with probability at least 1 - Gog^+0(klogN) ^ 2 . 7 ) 

for k G |0, KJ. The claim (12.7p is trivial for k = 0 since then pk = N and therefore (12.111 implies 
r*( 2 fc) ^ 1 deterministically. Now assume that (12.7p holds for some k G |0, Kj. Then Lemma [2.II 
applied with p = pk and M = 2 implies 

P(F*(zfc+i) ^4) ^ e-Gogi+0[kiogN) _ ^2.8) 

We may therefore apply Proposition 12.21 with z = z^+i and C = Clog? “ 0{klog N). Thus, we 
find that (|2.6p holds for z = Zk with probability at least 1 — Since \m\ ^ 1 

by (|1.15l) . we conclude that F*(zfc+i) ^ 2 with with probability at least 1 — 

This concludes the proof of the induction step, and hence of (12.6p for all Zk with k G |0, AT]. 

Finally, the argument may also be applied with ^ replaced by 2?, concluding the proof since 
g- 2 $ iog 2 $+ 0 (iog Ar)2 ^ g-^iog^ assumption. □ 

Remark 2.3. The induction in the proof of Theorem 1 1.1 1 is not a continuity (or bootstrapping) 
argument, as used e.g. in the works [MSI] on local laws of models with independent entries. The 
multiplicative steps p ^ p/2 that we make are far too large for a continuity argument to work, 
and we correspondingly obtain much weaker a priori estimates from the induction hypothesis. 
Thus, our proof relies on a priori control of F instead of the error parameters A^i and Aq used 
in [IMI]. The advantage, on the other hand, of the approach taken here is that we only have to 
perform an order log A^ steps, as opposed to the steps required in bootstrapping arguments. 
As evidenced by the proof of Theorem ll.il a logarithmic bound on the number of induction steps 
is crucial. An inductive approach was also taken in [13], where a local semicircle law without 
logarithmic corrections was proved for Wigner matrices with entries whose distributions are 
subgaussian. 
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It therefore only remains to prove Proposition 12.21 This is the subject of the remainder of 
the paper, which we now briefly outline. We follow the concentration/expectation approach, 
establishing concentration results on the entries of G (Section [3]) and computing the expectation 
of the diagonal entries (Section [5]). All of this is performed with respect to a conditional prob¬ 
ability measure, which is constructed for each fixed vertex. Roughly speaking, given a vertex, 
this conditional probability measure randomizes the neighbours of the vertex in an approxi¬ 
mately uniform fashion. It is model-dependent and has to be chosen with great care for all of 
the concentration/expectation arguments of Sections HHS to work. Its construction is easiest 
for the matching model, which we explain in Section [3l The constructions for the uniform and 
permutation models are given in Sections [6] and [7] respectively. 


3. Local resampling 

All models of random regular graphs that we consider are invariant under permutation of vertices. 
However, for our analysis, it is important to use a parametrization that distinguishes a fixed 
vertex. Without loss of generality, we assume this vertex to be 1. This parametrization has to 
satisfy a series of properties, which are given in Proposition 13.71 below. Using these properties, in 
Sections HHSl we complete the proof of Proposition 12.21 Loosely speaking, the parametrization 
allows us to resample the neighbours of 1, independently, and only changing a fixed number of 
edges in the remainder of the graph in a sufficiently random way. In this section, we describe 
the parametrization and prove Proposition 13.71 for the matching model. The parametrizations 
for the uniform and permutation models are discussed in Section [6] and [7] respectively. 

Random indices will play an important role throughout the paper. We consistently use the 
letters to denote deterministic indices, and x,y to denote random indices. 

3.1. Local switchings. Our basic strategy of local resampling involves randomizing the neigh¬ 
bours of the fixed vertex 1 by local changes of the graph, called switchings in the graph theory 
literature [35]. We use double switchings which involve three edges, as opposed to single switch¬ 
ings which only involve two edges. Both are illustrated in Figure [3T] 

Throughout the following, we use the following conventions to describe graphs. We consider 
general undirected graphs, which may have loops and multiple edges. We consistently identify a 
graph with its adjacency matrix A. The quantity Aij = Aji S N is the number of edges between 
i and j, and An S 2N is twice the number of loops at i. The degree of i is Aij, which will 
always be equal to d for all i. The graph A is simple if and only if it has no multiple edges or 
loops, i.e. Aij G {0,1} and An = 0 for all i,j. Sometimes we endow edges with a direction; we 
use the notation ij for the edge {i,j} directed from i to j. 

Let Aij denote the adjacency matrix of a graph containing only an edge between the vertices 
i and j, 

{Aij'jj^i — dn~5ji -|- 6iidji~. ( 3 . 1 ) 

To define switchings of a set of unoriented edges, it is convenient to assign directions to the 
edges to be switched. These directions determine which one of the possible switchings of the 
unoriented edges is chosen. We define the single switching of two edges rr, oo of A with the 
indicated directions to be the graph 

Trr_,aaj^At) '■= A -|- Aj-a -|- Aj-q — Aj-j. — A^a (3.2) 
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Figure 3.1. Solid lines show edges of a graph before switching, dashed lines edges after a single switching 
(left) and after a double switching (right). In our application, r is chosen to be 1, so that the switching 
connects the vertex 1 to a given vertex a. 


if {{''’tLi 0-1 Q,}\ = 4, and the graph Trr,aaiA) ■= A if Kr, r, a,a}| < 4. The double switching of the 
three edges rr, aa, bb of A with the indicated directions is defined to be the graph 


'^rr,aa,bbiA) — Tbr,aa{Trr,bbiA)') — A T ^ab T ^br ^rr ^bb ^aa (^A) 

if |{r, r, a, a, b,b}\ = 6 , and the graph Trr,aa. ,bbiA) ■■= A if \{r,r,a,a,b,b}\ < 6 . 

Our goal is to use switchings to connect the distinguished vertex 1 to essentially independent 
random vertices oi,..., that are approximately uniform in the sense of the next definition. 


Definition 3.1. A random variable x with values in [1, NJ is approximately uniform if the total 
variation distance of its distribution to the uniform distribution on is of order O (), 

i.e. ifEil^ix = *) - fI = 0{^)- 

To give an idea how approximately uniform random variables arise, consider a switching with 
r = 1 (to achieve our goal of connecting 1 to a given vertex a using a switching). For simple 
graphs, a necessary condition to apply the switching (13.3p is a 7 ^ 1. Choosing a uniformly with 
this constraint means that it is uniform on |2, A]. In particular, the total variation distance of 
its distribution to that of the uniform distribution is 

Throughout this paper, appears frequently as a bound on exceptional probabilities, and 
we tacitly use the estimates 

D ^ d, D ^ VdD < N, (3.4) 


which follow directly from (ll.5p - ()1.6p . as well as 

^ ^ $ 2 . (3.5) 

V d 

We use the following conventions for conditional statements. 

Definition 3.2. Let G be a a-algebra, B an event, and p G [0,1]. We say that, conditioned on 
G, the event B holds with probability at least p if¥(B\G) ^ p almost surely. Moreover, we say 
that, conditioned on G, the random variable x is approximately uniform if^j\^{x = i\G) “ fI ~ 
almost surely. 

The use of double switchings opposed to single switchings ensures that either condition (a) 
or (b) in the next lemma holds. These conditions will play an important role in Section |5l 
(That double switchings are in general more effective than single switchings is well known in the 
combinatorial context; see for instance [45] for a discussion.) 
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(3.6) 


Remark 3.3. Fix a d-regular graph A. For directed edges rl, aa, bb, dd, bb of A, we have 

~ '^rl,aa,bb{-^) — ^Id ~ ^ i 

where X is a sum of at most 8 terms zizAxy Explicitly, in the case 

\{l,r,a,a,b,b}\ = 1{1, r, a, a, 6,6}| = 6 (3.7) 

we have 

^ = ^ab ^bf ~ ^bb ~ ~ ~ ^bb + ^aa ■ ( 3 - 8 ) 

In particular, suppose that A is deterministic and the directed edges rl, aa, bb, dd, bb are random 
such that (13.7j) holds, a,a,d,d are approximately uniform, and, conditioned on a,a,d,d, the 
variables b, b, b, b are approximately uniform. Then for each term we have (a) the random 

variables x and y are both approximately uniform, or (b) conditioned on a, d, at least one of x 
and y is approximately uniform. 

We emphasize that when we say that x and y are approximately uniform, this is a statement 
about their individual distributions, and as such implies nothing about their joint distribution. 

The introduction of switchings that connect 1 to essentially independent random vertices 
ai,...,ad is simplest in the matching model, in which the different neighbours of any given 
vertex are independent, so that it suffices to consider a single neighbour of 1 at a time. In the 
next subsection, we explain in detail how this parametrization using switchings is defined for 
the matching model. 

We state the conclusion, Proposition 13.71 in great enough generality that it holds literally 
also for all of the other models, for which the more involved proofs are given in Sections [6HZ1 In 
the proof of Proposition 12.21 (given in Sections HHS]) , and therefore in the proof of Theorem ll.il 
we only use the conclusion contained in Proposition 13.71 and no other properties of the model. 
Hence, Proposition 13.71 summarizes everything about the random regular graphs that our proof 
requires. 

3.2. Matching model. The matching model was defined in Section ri.2l in terms of d independent 
uniform perfect matchings of |l,Ai]. We first consider one such uniform perfect matching, i.e. 
a uniform 1-regular graph. We denote by Sn the symmetric group of order N. For N even, 
denote by C Sn the set of perfect matchings of |1,A^], which (as explained in Section 
II. 2[) we identify with the subset of permutations whose cycles all have length 2; in particular 
TT = 7r“^ for TT G Mtv. For any perfect matching cr G M^v, we denote the corresponding symmetric 
permutation matrix by 

1 ^ 

(3-9) 

i=l 

Note that M(-) is one-to-one. 

Next, for i,j,k G |l,N'l, we define the switching operation : Mj^ Mjv through 

A4{Tijk{7r)) ■■= T^{i)i,jTr{j)M{k){M{7r)), (3.10) 

where we recall that r was defined in (13.3p . In particular, Tijk connects i to j (see Figure [3T]l 
except in the exceptional case \{i, j, k,7r{i),7r{j),iT{k)}\ < 6. 
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Lemma 3.4. Let vr he uniform over Mj\f, i G |1, A^J fixed, and a, b independent and uniform over 
[1, A^] \ {i}. Then Tiab{Tr) is uniform over M]\f, and 

{Tiab{'n-)){f) = a (3.11) 

provided that |{i, a, 6, 7r(i), 7r(a), 7r(6)}| =6. 

Proof. To prove that Tiab{T^) is uniform over Mjv, it suffices to check reversibility, i.e. that, for 
any fixed a, a' G M^r, 


= cr'lvr = cr) = ¥{Tiab{T^) = cr|vr = a'). (3.12) 

Given a, a' G Mjv, cr / o' , there is at most one pair {a,b) G ([1, A^]\{*})^ such that Tiab{o) = o', 
and such a pair exists if and only if there exists a (different) pair (a, b) such that Tiab{o') = o 
(see Figure ItT] (right) for an illustration). If no such pairs exist, both sides of (j3.1ip are zero. 
Otherwise, there exists precisely one pair (a, b) such that Tiab{o) = o', so that the left-hand side 
of (j3.11jl is equal to 1/(A^ — 1)^ because (a, 6 ) is uniformly distributed over {N — 1 )^ elements; 
the same argument shows that the right-hand side of (13.lip is also equal to 1/(A^ — 1)^, which 
concludes the proof of (j3.12p . Finally, (I3.11jl is immediate from the definition of Tiab- D 

The canonical realization of the probability space of the matching model is the product of d 
copies of the uniform measure on Mj\f. For our analysis, we instead employ the larger probability 
space LI := Lli X ■ ■ ■ X Lid where 


11^ := Mwx |2,iVlx |2,A^1, (3.13) 

also endowed with the uniform probability measure. Elements of Ll^ are written as (vr^, o^, b^). 
We set 0 = (vTi,.. .,TTd), Ufj, = {a^,h^), and 

:= ria^6^(vr^)- (3.14) 

By Lemma 13.41 ai,... ,cr^ are independent uniform perfect matchings of |1, A^J, and therefore 
the matching model is given by the adjacency matrix 

d 

T = Y.M{o^), (3.15) 

^l=l 

which is a random variable on the probability space LI. To sum up, rather than working directly 
with the probability measure on matrices that we are interested in, we use a measure-preserving 
lifting to a larger probability space, given by Ll^ Mjq —)> with (7r^,a^,6^) i-A- o^ := 

Throughout the following, we say that (ai,..., au) G |1, A^]'^ is an enumeration of the neigh¬ 
bours of 1 if 

d 

An = ^l(i = a^). (3.16) 

^i=i 

(Recall that, as explained in the beginning of Section [3l the vertex 1 is distinguished.) Defining 
Oifj. '■= o'^j(l)j we find that (ai,..., au) is an enumeration of the neighbours of 1. 
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3.3. General parametrization. Having described the probability space and the parametrization 
of the neighbours of 1 for the matching model, we now generalize this setup in order to admit 
other models of random regular graphs as well. 

Definition 3.5 (Parametrization of probability space). We work on a finite probability 
space 

n := 0 X C/i X • • • X c/d , (3.17) 

whose points we denote by (0, ui,..., Ud). Conditioned on 6 G Q, the variables ui,...,Ud are 
independent. For p € |l,d] we define a-algebras 

Ffj, := (3.18) 

Gfj, := a{e,ui,... ... ,Ud) ■ (3.19) 

We also define Fq ■= cr(0). 

In general, as in the case of the matching model in Section [3.21 the variable for p G [1, d] 
determines (with high probability given 0 G 0) the p-th. neighbour of 1. Note that we have 
introduced an artificial ordering of the neighbours of 1; this ordering will prove convenient 
in Sections HHS The interpretation of the cr-algebras (|3.18l) - (l3.19p is that determines all 
neighbours of 1 except the p-th one, and determines the first p neighbours of 1. 

Having constructed the probability space D, we augment it with independent copies of the 
random variables ui,..., Ud- 

Definition 3.6 (Augmented probability space). Let Lt be a probability space as in Defini¬ 
tion ro We augment Ll to a larger probability space D by adding independent copies of for 
each p G [l,d]. More precisely, we define 

D := 0 X C/i X • • • X C/d X C/i X • • • X C/d, (3.20) 

whose points we denote by {9,ui,... ,Ud,ui,... ,Ud). We require that, conditioned on 9, the 
variables ui,..., Ud, ui,... ,Ud are independent, and that and have the same distribution. 
On D we make use of the a-algebras defined by (j3.18p - (j3.19p . 

By definition, a random variable is a function X = X{9,ui,..., Ud, ui,..., Ud) on the aug¬ 
mented space Ll. Any function on D lifts trivially to an Fd-measurable random variable. Given a 
random variable X = X{9,ui,..., Ud, hi,..., Ud) and an index p G |1, d], we define the version 
X^ of X by exchanging the arguments and of X: 

Xf" := A(0,ui,...,u^_i,h^,u^+i,...,ttd,hi,...,h^_i,tt^,h^+i,...,'Ud). (3.21) 

Throughout the following, the underlying probability space is always the augmented space D. 
In particular, the vertex 1 is distinguished. However, since our final conclusions are measurable 
with respect to A = (Ajj), and the law of A is invariant under permutation of vertices, they also 
hold for 1 replaced with any other vertex; see in particular the proof of Lemma 15.41 below. 

Remark 13.31 and Lemma 13.41 imply the following key result for the matching model, which 
is the main result of this section. We state it in a sufficiently general form that holds for all 
graph models simultaneously; the proof for the other models is given in Sections [6HZ1 For the 
matching model, the parametrization in its statement and the corresponding random variables 
from (I3.22P were defined explicitly in Section [32} Ri below (13.131) . a* below (13.161) . and A in 

(ICT) . 
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Proposition 3.7. For any model of random d-regular graphs introduced in Section \1.B. there 
exists a parametrization satisfying Definition 1 ,9. ,51 augmented according to Definition \3.(A with 
Td-measurahle random variables 

tti,... 5 au^ ... 5 Old G |lj 5 (3.22) 

sueh that the following holds. 

(i) A is the adjacency matrix of the d-regular random graph model under consideration, and 
(ai,..., ad) is an enumeration of the neighbours of 1 in the sense of ()3.16n . 

(a) (Neighbours of 1.) Fix p. £ [l,d]. 

(1) Conditioned on Q^, the random variable a^ is approximately uniform. 

(2) Conditioned on Fq, with probability 1 — 0(-^) we have a^ = a^. 

(Hi) (Behaviour under resampling.) Fix fj, G |l,d]. 

(1) — A is the sum of a bounded number of terms of the form 2zAxy where x and y 

are random variables in |1, A^J. Conditioned on Q^, with probability 1 — the 

number of such terms is constant. Conditioned on Q^, for each term at least 

one of x and y is approximately uniform. 

(2) Conditioned on Fq, with probability 1 — 0 ('^ 7 ^) have 

A^-A = Aid, - Ala, + X , (3.23) 

where X is a sum of terms FAxy such that one of the following two conditions holds: 
(a) conditioned on Qy,, the random variables x and y are both approximately uniform; 
or (b) conditioned on Qy,ay,dy, at least one of x and y is approximately uniform. 
(Here we abbreviated Oy = Oy.) 

Proof of Proposition 1X71 matching model. The parametrization obeying Definition [331 
and the random variables (j3.22p for the matching model were defined in Section [321 We augment 
the probability space according to Definition 13.61 

The claim (i) follows immediately from Lemma 13.41 To show (ii) and (iii), we fix /i G |l,d], 
and drop the index p. from the notation and write for instance vr = vr^, a = Oy, and A = A^. 

First, we prove (ii). By definition, the random variable Oy is uniform on |2, A^J and hence 
approximately uniform on [1,A^], showing (ii)(l). By (j3.1ip . Uy = cr^(l) = Oy holds on the 
event |{1,7r(l), a, 7r(a), 6,7r(6)}| = 6 . The latter event has probability 1 — 0(;^) ^ 1 — 
conditioned on Qy, and hence in particular conditioned on Fq, which proves (ii)( 2 ). 

Next, we prove (iii). By the definitions (I3.14|) - (j3.15l) . 

A-A = M{T^-^{7r))-M{Tiai,{7r)). 

By the definition of T in (I3.10p and (13.3p , any application of T adds or removes at most 6 terms 
Axy, and therefore A — ^ is equal to a sum of at most 12 terms of the form FAxy, which proves 
the first claim of (iii)(l). 

To show the second claim of (iii)(l) and to show (iii)(2), we may assume that 

l{l,7r(l),a,7 r(a),6,7r(6)}| = 6 , |{1,7r(l), a, 7r(a), 6 ,7r(6)}| = 6 , (3.24) 
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since this event occurs with probability at least 1 — O(^) ^ 1 — 0(-^=) conditioned on (and 
hence also conditioned on J'o)- Under (j.S.24jl . we get 

A-A = M(Ti-^(7r)) - M{T,ab{7r)) = - T,i,,„,,b(M(7r)), (3.25) 

with r = 7r(l), a = 7r(a),6 = vr(6),a = Tr{a),b = 7r(6). As in Remark 3.3, we find that the 
right-hand side of (j3.25l) is 

Ala + — Aa7r(a) — (Ala + A^(i)^(ft) -|- Aft^(a) — — Aa7r(a)) ) 

from which the claim is obvious. □ 

3.4. Stability of the Green’s function under resampling. From now on we make use of the 
following notations for conditional expectations and conditional L^-norms. 

Definition 3.8. For any a-algebra Q, we denote by Eg = E( • \0) and Pg = P( • \0) the con¬ 
ditional expectation and probability with respect to Q. Moreover, we define the conditional L^- 
norms by 

WXhnG) ■■= (pe[i,oo)), 

ll-^llL°°(g) == sup{t > 0 : Pg(|A| > t) > O} . 

In particular, ||X||i;,p(g) is a G-measurable random variable, and 

Moreover, for any Td-measurable random variable X = X{6, ui,... ,Ud) we have 
\\^\\L<^{gA = raayi\X{e,ui,... . .Ud)\ ■ 

The following result is an important consequence of Proposition 13.71 for the Green’s function. 
It relies on the fundamental random control parameter 

= r^(z) := ||r(z)||i^oo(g^), (3.26) 

where we recall the definition of r(z) from ()2.4p . Also, we remind the reader that, according to 
Definition 13.61 a random variable (such as the index a; or y in the following lemma) is always 
dehned on the augmented probability space D, but the Green’s function is J-rf-measurable and 
does therefore not depend on ui,... , Ud- 

Lemma 3.9. Fix y e |l,dl. 

(i) For any i,j € |1, A] we have 

= Gij + 0{d-^/^r^T). (3.27) 

In particular, F^ = F -|- 0{d~^/'^T^T), and therefore F <C Vd implies F^ ^ 2r. 

(a) For random variables x,y such that, conditioned on and x, the random variable y is 
approximately uniform, 

= 0(r^$2)_ (3 28) 

An analogous statement holds with the roles of x and y exchanged, and with G replaced by 

G. 
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Assuming that F = 0(1), Lemma 13.91 (i) states that the Green’s function has a bounded 
differences property with respect to the u^: it only changes by the small amount = 

0(<I>) if a single is changed. Lemma 13.91 (ii) states that if one of its indices is random, then 
(conditioned on the L^-norm of the Green’s function is smaller (again by a factor $) than 
its L°°-norm. 


Proof. We start with (i). The resolvent identity (j2.2h implies 

= G,, + {d- ^ G,k{A - . (3.29) 

k,l 


By Proposition 13.71 (iii)(l), {A — A^)ki = 0 except for a bounded number of pairs (/c, 1), and the 
non-zero entries are bounded by an absolute constant. From this, we immediately get (|3.27l) . 

Next, we prove (ii). As in (I3.2nh . we may further augment the probability space to include 
another independent copy of which we denote by u^. From now on we drop the superscripts 
/i, and denote by X the version of an J'rf-measurable random variable X obtained by replacing 
with u^. On this augmented probability space, we introduce the fj-algebra := (t(^^,'S^). 
Then, since G is J^^-measurable (i.e. it does not depend on h^), we have Eg^/(G) = Eg /(G) 
for any function /. From (j3.27p . with replaced by we get 

G.J, = G^y + , 


and therefore 

\G^y\^ ^ 2\G^y\^+0{d-^V^y). 

Since, conditioned on Qy and x, the distribution of y has total variation distance G(-^^=) ^ 

O(-p) to the uniform distribution on |l,iV], and since |GxyP ^ ^ the Ward identity 

(|2.3h implies 


\Gxy\ ^ 2 


Im Gx 
Nr] 


+ 0{D-^Tl). 


Finally, by (13.271) . ImGa;^; ^ F^ + 0{D therefore 




xy\ 




2F^ 

Nr] 


+ 


Nr] 


+ oiD-^rt) = o{J^ + — + 




Nr] Nr] y/D 


which yields (I3.28p . 


□ 


4. Concentration 

In this section we establish concentration bounds for polynomials in the entries of G, with respect 
to the conditional expectation E_Fg. 

Proposition 4.1. Let z G C+ satisfy Ntj ^ 1 and let >0. Suppose that F = 0(1) with 
probability at least 1 — e“‘’. Then for any p = 0(1) and ... ,ipGp £ [1) A^l have 

Gnn • • • G^,J, - E^o [G,,,, • • • G,^,J = Om (4.1) 

with probability at least 1 — e“(Oog?)AC+ 0 (iogAf) _ 
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The rest of this section is devoted to the proof of Proposition 14.11 The main tool in its proof 
is the following general concentration result. 

Proposition 4.2. Let X be a complex-valued Xd-measurable random variable, and Yi,...,Yd 
nonnegative random variables such that Y^ is Q^-measurable. Let N satisfy d ^ Suppose 

that for all fi G |l,d| we have 

|X-Eg^X| ^ Y^, EgjX-Eg^X\^ ^ d-^Y^. (4.2) 

Suppose moreover that = 0(1) with probability at least 1 — e“^, and that Y^ ^ ^ 0 ( 1 ) ni^fYiost 
surely. Then 

X - = 0(0 , (4.3) 

with probability at least 1 — 

4.1. Proof of Proposition [121 To prove Proposition 14.21 we define the complex-valued martin¬ 
gale 

X^:=E^^X (//G[0,dl). (4.4) 

In particular, X^ = X and Xq = Ejr^X. By assumption, Y^ is bounded with probability 
least 1 — e“‘’. By the hrst inequality of (|4.2I) . we therefore get |^^+i — X^f = 0(1) with 
probability at least 1 — e”'’. If this bound held not only with high probability but almost surely, 
a standard application of Azuma’s inequality would show that X^ — Xq is concentrated on the 
scale Vd. This bound is not sufficient to prove Propositions 14. Ill4.21 which provide a significantly 
improved bound. Instead of Azuma’s inequality, we use Prokhorov’s arcsinh inequality, of which 
a martingale version is stated in the following lemma, taken from [291 Proposition 3.1]. Compared 
to Azuma’s inequality, it can take advantage of an improved bound on the conditional square 
function. 

Lemma 4.3 (Martingale arcsinh inequality). Let (Jyi)))^Q be a filtration of a-algebras and 
be a complex-valued -martingale. Suppose that there are deterministic constants 
M, So, si,..., Sd-i > 0 such that 


0^(L<d 


Then 


E{\Xd-XQ\^i) ^ 4exp - 


2y/2M 


arcsinh 


(s = 0,l,.. 

.,d-l). 

(4.5) 

&)) 

5 

(4.6) 


where S := •S/x- 

Proof. Since 


P(|Xd - Aol ^ ^ IP (l MXd - Xo)\ ^ lm{Xd -Xo)\^ 

it suffices to prove that any real-valued martingale X satisfying (14. obeys 


V2 


^{\Xd - Xo\ ^ f,) if 2exp -;—arcsinh — 




2M 


Mf 


2S 


(4.7) 


Hence, from now on, we assume that X is real-valued. 
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First, for all x G M, ^ 1 + x + x^(sinhx)/x and (sinhx)/x ^ (sinhy)/y if |x| ^ y. Using 
that (X^) is a martingale, it follows that for any A > 0, 

^ 1 + E^^(X^+ 1 -X^) 2 AsinhAM ^ l + s^AginhAM. 

Iterating this bound, using 1 + x ^ e^, it follows that 

^gA(x^-Xo) ^ exp (^Aginh(AM )5 

The estimate (j4.7p then follows by the exponential Chebyshev inequality with the choice A := 
^ arcsinh(-^), and an application of the same estimate with X replaced by —X. □ 

In order to exploit the fact that = 0(1) with high probability, we introduce a stopping 
time T. Let 7 ^ 1 be the implicit constant in the assumption of Proposition 14.21 such that ^ 7 
holds with probability at least 1 — e“‘’. We define 

r := min{// G |0,(i- 1] : ||y^+i^ 27 } , (4.8) 

and if the above set is empty we set r := d. By dehnition, r is an (J^^)-stopping time. The 
following result shows that r < d on an event of low probability. 

Lemma 4.4. Suppose that for all y G |l,d] we have P(L)i ^ 7 ) ^ and Y^ ^ ]\fO(i) Q^ijjiggi 
surely. Then 

P(T<d) ^ ^-C+OiiogN) _ 

Proof. For y G |0,d-l] set cj)^ := l(F),+i ^ 7 ). Then F),+i ^ 7 + and, by 

Minkowski’s inequality. 

Using a union bound, Markov’s inequality, logd = 0(log A^), 7 ^ 1 , and that = E0^ ^ 

e”'’ by assumption, we therefore get 

d—1 d—1 

P(r<d) ^ j;P(||y^+i||^2(^^)^27) ^ ^P(E^^iVO(i)^^^^ 2 ) 

/i=0 fi=0 

^ dN'^^^'>E4>^ ^ g-C+0(iogA)^ 

which concludes the proof. □ 

Since r is an (J^^)-stopping time, XJ^ •■= X^at is an (X^)-martingale. Because of Lemma 03] 
and using a union bound, it will be sufficient to study Xf instead of X^. The next result shows 
that Xf satisfies the assumptions of Lemma 14.31 

Lemma 4.5. For /x G |0, d — 1] we have 

|x;+i-x;| = 0(1), (4.9) 

e^jx;+i-x;| 2 = O(d-i). (4.10) 
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Proof. Set := l(r ^ /i + 1). Then is J-'^-measurable and 

x;+i -x; = MX,+1 - x^) = - E^^X). 

Note that, by definition, (/)^ = 1 implies that ||y^+i^ 27 = 0(1), and that, by indepen¬ 
dence, 

0^(X^+i - X^) = ct>^^T,UX - IEs,+i^) • (4.11) 

We now prove (14.91) . By the first bound of (|4.2I) . 

1-^ - ^g^+iX\ ^ y)j+i, 

and therefore 


4 >^l\X^l+l- x^\ ^ 0^Ejr^_^JX - Eg^^,^X| < ^ 27. 


In the last inequality, we used that = ())^Ejr^Y)t+i ^ ^ 27 = 0(1) 

since Y)t+i is ^^+i-measurable, by Holder’s inequality, and by the definition of (j)^. This completes 
the proof of (14.9|) . 

Next, we prove (|4.10l) in a similar fashion. By (14.lip . Jensen’s inequality for the conditional 
expectation and then using the second inequality of (|4.2p . we get 




/i+l 


-XI 


0^E^jX-Eg^^,X|2 ^ d-V^E^,y^'+i ^ 472^-1, 


as desired. □ 

Proof of Proposition 14.21 By Lemmas 0114.51 and ^arcsinh^ = Clog2^ + 0(1) for ^ > 0, 
we get 

P(|X - Ej-qXI ^ Ci) ^ P(|Xrf - Xq I ^ CCj + P{t <d) ^ Q-i^^og0A<:+O{iogN) ^ 
for a sufficiently large constant O. □ 

4.2. Proof of Proposition l4.ll Throughout the remainder of this section, we assume that Nr] ^ 1 
and 0^1. From Definitions 13.51(3^ we recall the cr-algebras and as well as the version 
X^ of a random variable X. In particular, we can express the conditional variance of an X^- 
measurable complex-valued random variable X as 

EgjX-Eg^Xp = iEgjX-X^|2. (4.12) 

The following result is the main ingredient in the verification of the second bound of (14.21) . 
For its statement, we recall the definition of F^ from (|3.26p . 

Lemma 4.6. We have 

EgJG.j - = 0{d-^Tl<^^) . (4.13) 

Proof. We abbreviate G = G^. Applying (I4.12D to G^-, we get 

Eg^lGjj — Eg^Gjjp ^ EgJ^Gij — Gij\^. (4.14) 
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Let X be the indicator function of the event of ^^-probability at least 1 — from Propo¬ 

sition O (iii)(l), and set x = 1 — x- Then the right-hand side of (I4.14h is bounded by 


- Gij)x\‘^ + \\Gij - Gij|||oo(g^)Eg^(x) • (4-15) 

To estimate both terms, we use that, by the resolvent identity and Proposition 13.71 fiiilflL 
there are a bounded (and possibly random) number £ of random variables (xi,yi),..., {x£,yi) 
such that 

N £ 

\Gij - Gij\ ^ {d- 1)-V2 ^ - A)kiGij\ = (d- 1)-G^ \GiccAJ • (4-16) 

k,l=l p=l 

We focus first on the second term of (j4.15p . By Proposition 13.71 liiilfll and ()4.16p . 

= 0{d-^/X) = 0(d-V2r3). 

By the definition of x and Proposition 13.71 (iii)(l), Eg^(x) = This implies 

that the second term in (I4.15P is bounded by the right-hand side of (|4.13l) . 

Next, we estimate the first term of (I4.15j) . By the definition of x and Proposition 13.71 (nil(1). 
the number £ in (I4.16P is constant on the support of x, and, conditioned on for each p G |1, ^|, 
at least one of Xp and yp is approximately uniform. Therefore 

1 ^ 

'^Oij.liGij — Gij)x\ ^ d — I I GjxpGy^jGjxqGy^j \, (4-17) 

p,q=l 

where, conditioned on Gp, for each (p, q), at least two of Xp, yp, Xq, yq are approximately uniform. 
We estimate two of the four factors of G or G by T^, including those without an approximately 
uniform index, and use the Cauchy-Schwarz inequality to decouple the remaining two factors of 
G or G, each of which has at least one approximately uniform index. Then using (|3.28p we find 
that each such term is bounded by 0(r®$^). Since the sum in (j4.17|) has a bounded number of 
terms, the claim follows. □ 

Proof of Proposition [47T1 We verify the assumptions of Proposition 14.21 Given p = 0(1), 
set Yp, ■■= for a sufficiently large constant Gp. By definition, Yp is ^^-measurable. 

Moreover, by assumption, P = 0(1) with probability at least 1 — e“^. Hence, Lemma 13.91 
(i) implies that P^ ^ 2r = 0(1) with probability at least 1 — e“^, so that Yp = 0(1) 
with probability at least 1 — e”'’. Moreover, the trivial bound (El]) and Np ^ 1 imply 
Yp = 0{p~‘^~P) = 0{N^^P) = We conclude that Yp satisfies the conditions from the 

statement of Proposition 14.21 

We first complete the proof for p = 1. Let X := Then, by ()3.27|1 and $“3 ^ d^G^ 

|X-Eg^X| < 4.-i|G,,-Eg^G,,| = ^-^0{d-^Grl) = 0(r2) ^ Yp, (4.18) 

assuming that the constant Gp was chosen sufficiently large. This establishes the first estimate of 
(14.21) . The second estimate of (14.2p follows from Lemma 14.61 Therefore Proposition 14.11 follows 
from Proposition 14.21 
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Next, we deal with the case of general p. For k G |l,p] abbreviate Qk '■= Gi^j^ and consider 
X := ^~^Qi ■ ■ ■ Qp. By telescoping, <1>(X — Eg^X) is equal to 

P r . 

\Qi ■ ■ ■ Qk-l{Qk — '^g^Qky^g^{Qk+l • • • Qp) — Ql • • • Qk-l^g^{{Qk — ^g^Qk){Qk+l • • • Qp)) ; 

k=l 

and therefore 


\X-Eg^X\ ^ <l>-^TP-^^[\Qk-Eg^Qk\+EgJQk-Eg^Qk\ . 

k=l 

Using (13.2711 ■ we therefore conclude that \X — Eg^X\ ^ (after choosing Cp large enough). 
Moreover, since (oi + • • • + a 2 p)‘^ ^ (2p)^(a^ + • • • + a|p), by the conditional Jensen inequality 
and Lemma 14.61 we find 

Eg^\X-Eg^X\^ ^ 0(p2)cl,-2r2p-2maxEgJG,,-Ee^G,,|2 < 

which is bounded by dr^Y^ (after choosing Cp large enough). The claim now follows from 
Proposition 14.21 □ 


5. Expectation 

In this section we prove Proposition 12.21 We use the spectral parameters 

zq = E + \r]Q, z = E + ir], ^‘^/N < r/o < r? < IV. (5.1) 

Fix zq as in (I5.1jl . To prove Proposition 12.21 we assume that D ^ and that 

P(F*(zo)^ 7 ) = P(max{F(£'+ ir/) : 77 ^ ? 7 o} ^ 7 ) ^ e“'’ (5.2) 

for some constant 7 = 0{1). Recall the function E = Ez from (ll.ll|) and $ from (jl.lOll . To 
prove Proposition 12.21 it suffices to show that, with probability at least 1 — 

max|Gjj —m| = 0(T(^4>)), (5-3) 

i 

max|Gjj| = 0((^<h), (5-4) 

i¥=j 


for any z satisfying dSH). 

The proof of (j5.3l) - (j5.4l) proceeds in the following steps: 

(i) Estimate of s — m, where s is the Stieltjes transform (ll.lbjl of the empirical spectral 
measure, and m the Stieltjes transform (jl.9p of the semicircle law. 

(ii) Estimate of Gu — m . 

(iii) Estimate of Gij for i j. 

Step (i) represents most of the work. Throughout this section we make the assumption (15.2j) . 
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5.1. High probability a priori bounds. For the proof of Proposition 12.21 we use the following 
convenient notion of high probability. 


Definition 5.1. Given a parameter t > 0, an event H holds with t-high probability, abbreviated 
t-HP, i/P(H'=) ^ e-t+0{\ogN)_ 


In the nontrivial case t S> log N , the notion of t-high probability is stronger than the stan¬ 
dard notion of high probability (and in fact implies what is occasionally called overwhelming 
probability). By definition and a union bound, an intersection of many events that each 

hold with t-high probability holds with t-high probability. Moreover, if H holds with t-HP then 
]Ej-q1(H^) ^ N~^ with t-HP for any constant k > 0. Indeed, by Markov’s inequality, 

P(Ej-ol("^) > 1/A^^) ^ ^ Q-t+0[iogN) _ 

From now on, these properties will be used tacitly. 

Furthermore, from now on, the parameter t in Definition 15.11 will always be 

t := (^log^AC (5.6) 

with () and ^ the parameters given in the assumption of Proposition 12.21 Then, for any z as in 
(EH), we get from the assumption (j5.2p and Proposition 14. II that, with t-HP, for all deterministic 
i,j,k,l,m,n G [l,iV|, 

\G^j\ = 0(1), = Gi, + 0{(^), (5.7) 

and 

^ToiGijGkl) = GijGkl + 0{^^) , Ejr^{GijGklGmn) = GijGklGmn + 0{^^) . (5.8) 

To prove Proposition 12.21 we need to show that (|5.3p - ()5.4l) then also hold with t-HP. 

5.2. Derivation of self-consistent equation. In this subsection we derive the self-consistent equa¬ 
tion, (I5.35P below, which will allow us to obtain estimates on the entries of G and hence prove 
Proposition 12.21 The following lemma is, in combination with the concentration bounds (j5.7p - 
(15.81) . the main estimate in its derivation. For its statement, recall from Proposition 13.71 that 
(oi,..., ad) is an enumeration of the neighbours of 1. For the following we introduce the abbre¬ 
viation 

( 5 . 9 ) 

i 

so that, under i is regarded as a uniform random variable that is independent of all other 
randomness. With this notation, we may express the the Stieljes transform (I1.16P of the empirical 
spectral measure as s = eWG**. 

Lemma 5.2. Fix fi G [l,d]. Given z G C+ with Nrj ^ 1, suppose that F = 0(1) with t-HP. 
Then for all fixed j,k,l G [1, iV], 

ETo{Gc,,j-E^^Gij + {d-l)-^/hGij^ = 0(d-i/2^), (5.10) 

E^,(Gki[Ga,j -E^^Gij + {d-l)-G^sGij)) = 0(d-i/2<^), (5.11) 


with t-HP. 
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Recall from (j3.2ip that is the version of a random variable X = X{9, ui, ..., Ud, ui,..., Ud) 
obtained from X by exchanging its arguments and u^. Throughout this section we make use 
of the indicator function 

X^i •= where := li«,, = a,,) 1 (claim (13.231) from Proposition 13.71 holds) . (5.12) 

Note that Xij. = X^- Moreover, by Proposition 13.71 liiR2) and (hi) (2) as well as a union bound, 
we have Ej-o(x^) = 1 - O(^). 

For brevity, given a fixed index ^ G |l,(i], we often drop sub- and superscripts and write 
simply 

a = a = a = , y = , A = , G = &. (5.13) 

(As in Proposition 13.71 we always abbreviate = d^.) 

The following lemma provides several elementary bounds on the Green’s function. It is the 
main computational tool in the proof of Lemma 15.21 

Lemma 5.3. Given z G C+ with Nr] ^ 1, suppose that T = 0(1) holds with t-HP. Fix // G [1, d}, 
and use the abbreviations (j5.13p . Then the following estimates hold with t-HP. 


(i) For all j G |1, we have 


E^„(G,,) = E^„eW(G,,) + 0 (^), 

(5.14) 

^MGaaGjj) = E^,E^^{GuGjj) + 0{^). 

(5.15) 

(a) For all i,j,k,l,m,n G {1, NJ we have 


E^,{xG,j) = E^,{G,j) + 0 {^), 

(5.16) 

E^.ixGijGki) = E^,{GijGki) + 0 {^), 

(5.17) 

Ej^QixGijGklGmn) — Ejr^^{GijGklGmn) + . 

(5.18) 

Analogous statements hold if some factors G are replaced with G. 


(Hi) For any i,j,k,l,m,n G {1, NJ we have 


Ej^oiGijGki) = Ejr^{GijGki) + 0{^) , 

(5.19) 

E^,{GijGklGmn) = E^,{GijGklGmn)+0 {^). 

(5.20) 

(iv) If (a) conditioned on and d, the random variable x is approximately uniform, or (b) 

conditioned on Q^, the random variable y is approximately uniform, then 

E^,{Gda,Gyi) = 0(d>), 

(5.21) 

E^,{G,jGd,,Gyi) = 0(^>). 

(5.22) 

Proof. Fix u G [1. d]. and. as in the statement of the lemma, use the shorthand notation (15.131). 
Denote by (f the indicator function of the event F^ ^ 27 , and set (f = 1 — <j). By definition, cj) 
is Gy,-measurable. Bv (j3.27l). |F^, < 27 } C |F ^ 7 ). so that, bv assumption. </> = 1 with t-HP. 
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for a constant 7 = 0(1). In particular, as noted around (I5.5p . for any constant /c, the event 

^ ^ 7 ^} 

(i) We show (I5.14p : the proof of pS.lSp is analogous. Since, conditioned on Jq) h and G are 
independent, and the total variation distance between the distribution of a and the uniform 
distribution on |1, A^] is 0(-^7=), 

E^,{Gaj) = E^.iGajcP) + 

= + O(^) = E^3EW(G,,) + O(^) , (5.23) 

with t-HP. 

(ii) We show (I5.17P : the proofs of p5.16p and p5.18p are analogous. Since x ^ 1) 

Ej-o((l - x)GijGki) = E^„(<^(1 - x)GijGki) + 0{E^,{r ,-^)). (5.24) 

The first term is bounded by 0(E_7,(1 — x)) = by the definition of x and Propo¬ 

sition [3]71 The second term is also O () with t-HP, as observed at the beginning of the 
proof. 

(iii) We show p5.19p : the proof of (j5.20l) is analogous. Since r]~‘^Ejrg^ = O(^) with t-HP and 
4>\Gij\ = 0(1), we get from (I3.27P that 

EjToiGijGki) = Ej-g{(f>GijGki) + 0{^) = Ejr^{(f>GijGki) + 0{^) 

= E^„(G,,Gfc0+O(^), (5.25) 


with t-HP. 

(iv) We show (15.2111 : the proof of (j5.22p is analogous. Under assumption (a), the Cauchy-Schwarz 
inequality and (I3.28P imply 

\Ej^Q{GaxGyl)\ ^ Ejr^(j)\GaxGyl\ + r]~‘^Ejr^^ 

^ 0{E^,4>EgJGax\^)G^ + = 0(d>), (5.26) 

with t-HP, where we used (j)EgJGyi\‘^ ^ (/>r^ ^ 47^ = 0(1). Similarly, under assumption (b), 

\EMGaxGyi)\ ^ O(E^Q0Egj4i|2)V2 + ^-2]E^^^ = 0(<h), (5.27) 

with t-HP, where we again used (15.5p . This completes the proof. □ 

Proof of Lemma 15.21 The proofs of both estimates are analogous, and we only prove (j5.10p . 
Throughout the proof, we use Lemma [5^ repeatedly, and estimate ^ Since /r G |l,d] 

is hxed, we also use the abbreviations (I5.13P in the remainder of the proof, and use the indicator 
function x = X/^ dehned in (I5.12p . By definition, conditioned on J-q, the random variables u and 
u are identically distributed, so that Gaj and Gaj are also identically distributed. (Recall the 
definition p3.21l) and the convention (j5.13p .i Thus, by (I5.16p . and since a = a on the support of 
X (by definition (I5.12jl of x), we obtain 

E^,{Gaj) = E^.ixGaj) + 0(d-i/2$) = E^,ixG~aj) + 0(d-i/2^), (5.28) 
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with t-HP, where in the last step we also used that X = X^- By (I5.14p and (I5.16p . with t-HP, 


Ej-oEW(G,,-) 

= E^ffxGaj) + 0{d-^/H). 

(5.29) 

This implies, with t-HP, 

E^ffGaj-E^^Gij) 

= E^,x{Gaj - Gaj) + 0{d-^/^<I>) . 

(5.30) 


By the resolvent identity, G — G = {d — 1) ^/‘^G{A — A)G, and therefore by Proposition 13.71 
(iii)(2), on the event {x = 1} we have 

G-aj - Gaj = (d- l)-V2(_Gs-(5y + s) , (5.31) 

where S a sum of a bounded number of terms of the form zizGaxGyj with random variables 
X and y such that at least one of the following two conditions is satisfied: conditioned on 
and d, the random variable x is approximately uniform, or, conditioned on Gy, the random 
variable y is approximately uniform. (For example, S contains the term GaaGij corresponding 
to {x,y) = (a, 1). Conditioned on Gy, the random variable x = a is approximately uniform and 
independent of a, so that x is approximately uniform conditioned on Gy and d.) Therefore, by 
(j5.17p . (I5.19p . and (I5.2ip . we get 

\E^,{xS)\ ^ E^„|5| + 0(ci)) = 0(cl.), (5.32) 

with t-HP. Similarly, by (15.171) . (15.191) . and (I5.15p . we get 
E^.ixG-aaGij) = E^,{GaaGy) + 0{<^) 

= E^,{GaaGy) + 0{<^) = E^,E^^{GiiGij) + 0{<^), (5.33) 


with t-HP. From ()5.3np - (j5.33p . we conclude that 

(G„,- - E^^Gi,) = {d- l)-i/2 (-E^„(EWGiiGi,) + (9($)) , (5.34) 

with t-HP. Since EWGi* = s, we obtain (|5.10p . The proof of (|5.1ip is analogous, using (I5.18P 
instead of (j5.17p . (j5.20p instead of ()5.19p . and (j5.22p instead of ()5.2ip . □ 

The main idea of the proof of Lemma [5.2l is (I5.30p : the left-hand side is a difference of Green’s 
functions with different indices, while the right-hand side is (up to a small error) a difference of 
Green’s functions with the same indices but the first Green’s function is computed in terms of 
a switched graph. 

We now have all of the ingredients to derive the self-consistent equation for the diagonal 
entries of G. 

Lemma 5.4. Given z G C+ with Nrj ^ 1, suppose that (15.7p - (15.81) hold with t-HP. Then, we 
have with t-HP, for all j G |l,iV], 

1 -|- (s -|- z)Gjj = 0((1 -|- |^|)c^4*). (5.35) 

In particular, with t-HP, 

l + sz + s"^ = (9((1-h l^DC^). (5.36) 
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Proof. The event that (|5.35l) holds is measurable with respect to A. By invariance of the law 
of A under permutation of vertices, and a union bound, it therefore suffices to establish (j5.35p 
for j = 1 only. Then (|5.36p follows by averaging (|5.35p over j. 

To show (j5.35p with j = 1, we make use of the larger probability space hi from Definitions l3.5l - 
13.61 where the vertex 1 is distinguished. By (I5.7p - (l5.8p . it is sufficient to show that, with t-HP, 

l + zEjr^Gii = —IEj-q (sGii) + 0(<1>). (5.37) 

To show (j5.37p . we use that by {H — z)G = I and (II.7p . with (j5.9p . 

1 + zGii = EE( 

i i /i=l ^ ^ 

d 

= (5.38) 

^=1 

Taking the conditional expectation on both sides of (I5.38P and using Lemma [521 we get 

l + zE^,Gn = -^E^„(sGii) + 0($) (5.39) 

with t-HP. This implies (I5.37P and therefore completes the proof. □ 

Under the assumptions of Proposition 12.21 the statement of Lemma (5.41 mav be strengthened 
as follows. 

Lemma 5.5. Let zq be as in dST]) and suppose that (I52D holds. Then with t-HP the estimates 
(I5.35p - (j5.36l) hold simultaneously for all z as in (15.ip . 

Proof. Set 

m:=po + l/N\ le[0,iV5], (5.40) 

and Zi := E + ipi. Since (j5.7l) - (l5.8p hold uniformly with t-HP for any rj ^ r/o, by Lemma [5.4l and 
a union bound, (j5.35p - (l5.36p hold simultaneously at all zi with I G |0,A^^], with t-HP. Since 
is a l/iV^-net of [? 70 ) Vo + Gij is Lipschitz continuous with constant 1/?]^ ^ N'^, the 

claim follows. □ 

5.3. Stability of the self-consistent eqnation. In Lemma [53] we showed that, with t-HP, 

s^ + sz + l = 0((l + |z|)C4>). (5.41) 

It may be easily checked that the Stieltjes transform of the semicircle law ()1.9p is the unique 
solution m : C+ —)• C+ of the equation 

+ mz -|- 1 = 0 . (5-42) 

To show that m and s are close, we use the stability of the equation (15.4211 . in the form provided 

by the following deterministic lemma. The stability of the solutions of the equation (15.421) is a 
standard tool in the proofs of local semicircle laws for Wigner matrices; see e.g. |21] . Our version 
given below has weaker assumptions than previously used stability estimates; in particular, we 
do not (and cannot) assume an upper bound on the spectral parameter E. 
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Lemma 5.6. Let s : C+ —)• C+ be continuous, and set 


R := + sz + 1. (5.43) 

For E G W, r]Q > 0 and rji ^ 3 V r]Q, suppose that there is a nonincreasing continuous function 
f [0,1] such that \R{E + i? 7 )| ^ (1 + + ir]\)r{rj) for all g G [ 770 , 771 ]. Then for all 

z = E + \g with 77 G [t?o, ??i] we have 


s — m\ = 0{F{r )), 


(5.44) 


where F was defined in dm]). 


Proof. Denote by m and m the two solutions of (I5.42|) with positive and negative imaginary 
parts, respectively: 

+ ... -Z-VW^^ 


m = 


m = 


2 2 

where the square root is chosen so that ImTn > 0. Hence, Im z"^ — 4 > 77 and consequently 
ImTu < — 77 ; we shall use this bound below. Note that m and m are continuous. Set v '■= s — m 
and V := s — m. Since 

—z ± — 4 + 4ii 


s = 


and since for any complex square root and tc, C G C we have 

I a/tu + c - \Ac| A I a/tt; + C + Vwl ^ 


we deduce that 


k; A 7 ) < 


m\ 


Vh 


A^/\R\ ^ 3F(r). 


In the last inequality, we used that |ii| ^ (l + |z|)r and that, for any r G [0,1], 

AI + NK < 3F(r). 

Vk^-4| 


(5.46) 


(5.47) 


(5.48) 


(5.49) 


The proof is divided into three cases. First, consider the case (1 + |z|)r( 77 ) ^ \z‘^ — 4|/16. 
Then, using (15.491) and the fact that \v — v\ = ■\/\z^ — 4|, we get 


|7| ^ [u] + \/\z^ — 4| = juj + — AA/(l + |z|)r| = |7;| + 0(F(r)), (5.50) 

Vvk^-4| ) 

and hence juj ^ juj A IDj + 0{F{r)) = 0{F{r)) by ()5.48p . 

Next, consider the case 77 ^ 3. Then, on the one hand, by (I5.48|) and the assumption r G [0,1], 
we have juj A |{)| ^ 3F{r) ^ 3. On the other hand, since Ims > 0 and ImTu, < —77 ^ —3, 


lul ^ jlms —ImTuj > jlmThj > 3, 

and together we conclude that | 7 ;| > juj, so that the claim follows from (I5.48p . 


(5.51) 


29 


































Finally, consider the case (1 +|2:|)r(77) < —4|/16 and r] < 3. Without loss of generality, we 

set r]Q = rj. Since r is nonincreasing and \z^ — 4| increasing in r/, and since (1 + \z\) ^ 4(1 + |zo|) 
for T] £ [r/o, 3], we then have (1 + \z\)r{r]) < \z'^ — 4|/4 for all t] £ [tjq, 3]. Therefore 

\v-v\ = \/\z^ - 4| ^ t A ^/(l + |z|)r [ for all r/£ [r?o,3] . (5.52) 

Vvk -4| / 

On the other hand, by (15.481) . and since |i?| ^ (1 + \z\)r, 

l^l A |{)| ^ ^ \/(l + \z\)r for all r] £ [r]o, 3] . (5.53) 

VP -4| 

By continuity and (j5.52^ - 115. 53 j) . it suffices to show |u( 2 ;)| < |i)(z)| for some r/ £ [%) 3], since then 
\v{z)\ < \v{z)\ for all rj £ [%)3]. Since we have already shown |u( 2 ;)| < |'0(2:)| for r] = 3, the proof 
is complete. □ 

5.4. Proof of Proposition 12.21 We now have all the ingredients we need to complete the proof 
of (|5.3l) - (l5.4p under the assumption (15.2p . and hence the proof of Proposition 12.21 

Proof of (j5.3p . Let zq = E + irjo £ D he given, where D was defined in p2.5p . Set r/i = N. 
Lemma 15.51 shows that, with t-HP, for all 77 £ [ 770 , 771 ], the function s satishes (|5.43l) with 

\R{z)\ V (l + |2;|)r(77), (5.54) 

where r{r]) ;= + ir]) and C is some large enough absolute constant. Hence, r is decreasing 

in 77 and, since <C 1 by assumption, we have r £ [0,1]. From Lemma 15.61 it therefore follows 
that |m — s| = 0(F(C'^4>)) = 0(F(^$)) for all 77 £ [ 770 , 771 ], with t-HP. 

Having determined s, we now estimate Gjj — m. By (I5.35P and since s = m + 0(F(^4>)), we 
find 

l + {z + m)Gjj = 0(F(e$))G^y+ 0((l + |z|).e^) (5.55) 

with t-HP. From (jl.9l) it is easy to deduce that z + m = —1/m and (1 + [ 2 :|)|m| = 0(1). Hence, 
by (I5.55P and since Gjj = 0(1) with t-HP, 

777 - 0 ,, = o{Em)G,,+om = o{Fm) ( 5 . 56 ) 

with t-HP, as claimed. □ 

The off-diagonal entries of the Green’s function can be estimated using a similar argument. 
Proof of (15. 4p . From {HG)ij = {{H — z)G)ij + zGij = 5ij -|- zGij, it follows that 

Oi 2 (HG)ii-Gii(HO)i 2 = Gi 2 , (5.57) 

and therefore 

Gi 2 = Gi 2 Hi,G il - On ^ HiiG ,2 . (5.58) 

i i 

As in (I5.38p . using p5.58l) and (|1.7I) . we hnd 

d 

Gi2 = -(d-l)-'/'j;EW(Gii(G«^2-0,2)-Gi 2 (G„^i-Gii)). (5.59) 

it=i 
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Figure 6.1. Each of the three diagrams illustrates the subgraph incident to a set S C E of three distinct 
edges (IS'! =3). We draw the edges of S with solid lines and the edges of E \ S' with dotted lines. In the 
left and centre diagrams we have I{E, S) = 0: in the left diagram |[S]| < 6, while in the centre diagram 
|[S]| =6 but i?|[s] ^ S. In the right diagram we have I{E, S) = I. 


Using (jS.lip we therefore get, with t-HP, 

d 

^ToGi 2 = -id - 1)-^ (^115012 - GusGu) + 0(^>) = 0(^>) • (5.60) 

By Proposition 14.11 therefore G 12 = E^pGi 2 + = 0(^4)) with f-HP. Again, by symmetry 

and a union bound, the claim then holds uniformly with 12 replaced by ij. □ 

Summarizing, we have proved that, assuming (15.2p . the estimates (|5.3p and (15.41) hold with 
t-HP. Hence the proof of Proposition 12.21 (and consequently of Theorem II.ip is complete. 

This proof of Theorem 11.11 relies on Proposition [321 which we proved for the matching model 
in Section 13.21 In order to establish Theorem o for the uniform and permutation models, we 
still have to prove Proposition 13.71 for these models. This is done in Sections [6] and [71 which 
constitute the rest of the paper. 


6. Uniform model 

In this section we prove Proposition 13.71 for the uniform model. We identify a simple graph on 
the vertices [1, A"] with its set of edges E, where an edge e £ E is a subset of |1, A] with two 
elements. The adjacency matrix of a set of edges E is by definition 

M(A) := ^ Aij, (6.1) 

{*d}6E 


where Ajj was defined in dST]). Note that M(-) is one-to-one, i.e. the matrix M{E) uniquely 
determines the set of edges E. For a subset 5 C A of edges we denote by [5] := UeeS ^ 
of vertices incident to any edge in S. Moreover, for a subset B C [1, A] of vertices, we define 
E\b ■= {e £ E : e d B} to be the subgraph of E induced on B. 

6.1. Switchings. For a subset S d E with \S\ = 3 we define the indicator function 
/(A, 5) := 1(A|[5] is 1-regular) = 1(|[S’]| = 6, A|[ 5 ] = S'). 
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The interpretation of I{E, S') = 1 is that S is a switchable subset of E, i.e. any double switching 
of the three edges S results again in a simple d-regular graph; see Figure IHTTl A switching of the 
edges S may be identified with a perfect matching of the vertices [S]. There are eight perfect 
matchings S' of [S] such that S n S' = 0. We enumerate these matchings in an arbitrary way as 
S' with s G [1, 8 ], and set 

Ts,s{E) := (F;\S)US', ( 6 . 2 ) 

and say that Ts^siE) is a switching of E. (Compare this with Figure [3T] (right) in which one 
such perfect matching is illustrated.) Note that there are rr, aa, bb depending on (S, s) such that 

M{Ts,s{E)) = Trr,aa,bb{M {E)) , (6.3) 

with the right-hand side defined by (13.3p . This correspondence will be made explicit later. The 
definition (j6.2l) implies, for I{E,S) = 1, that Ts^s{E) is a simple d-regular graph, and that 

E\S = TsAE) n{E\S). (6.4) 

Next, take two disjoint subsets Si, S 2 C E satisfying I{E, Si) = I{E, S 2 ) = 1 and [Si]n[S 2 ] = 
{!}. Thus, we require the sets Si and S 2 to be incident to exactly one common vertex, which 
we set to be 1; see Figure E21 Then S 2 <Z E \ Si and Si <Z E \ S 2 , and (16.4p implies that the 
two compositions Tsj^^si{Ts 2 ,s 2 {E)) and Ts 2 ,s 2 {Esi,si{E)) are well-defined and coincide, 


TSi,sATs2,S2{E)) = Ts2,S2{TSi,Si{E)) . (6.5) 

Let S satisfy I{E,S) = 1 and 1 G [S]. The map Ts,s{E) switches the unique edge {l,i} G S 
incident to 1 to a new edge {1, j} ^ S with j G [S]. Our next goal is to extend this switching to 
a simultaneous switching of all neighbours of 1. As already seen in (16.51) . simultaneous multiple 
switchings are not always possible, and our construction will in fact only switch those neigh¬ 
bours of 1 that can be switched without disrupting any other neighbours of 1. The remaining 
neighbours will be left unchanged. Ultimately, this construction will be effective because the 
number of neighbours of 1 that cannot be switched will be small with high probability. 

Let {ei{E),... , ed{E)) be an enumeration of the edges in E incident to 1, and denote by 

SAE) := {5cU:e^(F;)G5,|5|=3,l^eforeGS\{e^(F;)}} ( 6 . 6 ) 

the set of unordered triples of distinct edges in E containing e^E) and no other edge incident 
to 1. Conditioned on E, we define a random variable (S, s), where S = (5i,..., Sd) and s = 
(si,..., Sd), uniformly distributed over Si{E) x • • • x Sd{E) x |1,8]'^. In particular, conditioned 
on E, the random variables (^i, si),..., {Sd, Sd) are independent. 

For G |l,d] we define the indicator functions 

A = IAE,S) ■.= I{E,SA, (6.7) 

Jfi = JA^) •= l([5';,] n [S'i.] = { 1 } for all 1 ///r) , ( 6 . 8 ) 

and the set 

W = W{E, S) := {// G 11, dj : ^(U, S) J^(S) = 1} . (6.9) 

Their interpretation is as follows. On the event {/^ = 1}, the edges 5^ are switchable in the 
sense that any switching of them results in a simple d-regular graph. The interpretation of 
{Jf 2 = 1 } is that the edges of S^ do not interfere with the edges of any other S^, and hence any 
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Figure 6.2. Each of the two diagrams illustrates the subgraph incident to Si and S 2 for two sets of 
three edges S*!, S '2 C E satisfying I{E, Si) = I{E, S 2 ) = 1. The edges of Si U S 2 are drawn with solid 
lines and the edges of E\ {Si U S' 2 ) with dashed lines. The vertex 1 is drawn using a white circle. In 
terms of the indicator functions defined in (I6.8|l . in the left diagram we have Ji{E, S) = J 2 {E, S) = 0 
since [^i] fl [5'2] ^ {!}, and in the right diagram Ji{E, S) = J 2 {E, S) = 1. 


switching of them will not influence or be influenced by the switching of another triple of edges: 
on the event {J^ = 1}, (16.Sp implies that commutes with for all u / /r. The set W 

lists the neighbours of 1 that can be switched simultaneously; see Figure [621 
Let /ii,..., /ifc, where /c ^ d, be an arbitrary enumeration of IT, and set 

Ts,s{E) ■■= ( 6 . 10 ) 

By (16.5p . the right-hand side is well-defined and independent of the order of the applications of 
the Ts^^Sfj,- Equivalently, in terms of adjacency matrices, Ts^s{E) is given by 

M{Ts,s{E))-M{E) = Y,{M{Ts,,s,{E))-M{E)), (6.11) 

fieW 

where we used that, by construction of W, any switchings ^ with £W do not interfere 
with each other, so that M{Ts^,s^Ts,,sA^)) “ M{Ts,,s,{E)) = M{Ts^,s^{E)) - M{E). 

The following result ensures that the simultaneous switching leaves the uniform distribution 
on simple d-regular graphs invariant. For this property to hold, it is crucial that, as in (16.6p . we 
admit conhgurations S that may have edges that cannot be switched. The more naive approach 
of only averaging over configurations S in which all edges can be switched simultaneously does 
not leave the uniform measure invariant. The price of admitting conhgurations S that do not 
switch some neighbours of 1 is mitigated by the fact that such conhgurations are exceptional and 
occur with small probability, i.e. conditioned on E, /^(E, S) = T/i(S) = 1 with high probability. 

Lemma 6.1. If E is a uniform random simple d-regular graph, and (S, s) is uniform over Si{E) x 
• • • X Sd{E) X |1,8]'^, then Ts^s{E) is a uniform random simple d-regular graph. 

Proof. It suffices to show reversibility of the transformation E 1 — )• Ts^s{E) with respect to the 
uniform measure, i.e. that for any hxed simple d-regular graphs Ei,E 2 we have 

F{Ts,s{E) = E 2 \E = El) = P(rs,s(E) = Ei\E = E 2 ) . (6.12) 

Note that (S,s) is uniformly distributed over iSi(£'i) x • • • x iSrf(Fii) x |1,8]'^ on the left-hand 
side and over Si{E 2 ) x • • • x Sd{E 2 ) x |1, 8]'^ on the right-hand side. 
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Figure 6.3. The lines (solid and dashed) illustrate the edges of E^. The solid lines are the edges of 
El and the dashed line the edges of E 2 . The white circle represents the vertex 1. Note that E/\\^ 2 ,...,n} 
consists of |hF| = 2 disconnected subgraphs, each of which has 5 vertices and 4 edges. These subgraphs 
are given by the encircled regions with the vertex 1 and its incident edges removed. The sets are 
given by the encircled vertices including 1. 


First, given two (simple d-regular) graphs Ei,E 2 , we say that E 2 is a switching of Ei if there 
exist (S,s) such that E 2 = Ts^s{Ei), and note that Ei is a switching of E 2 if and only if E 2 
is a switching of Ei. If these conditions do not hold, then both sides of (I6.12p are zero. We 
conclude that it suffices to show (I6.12|) for the case that E 2 is a switching of Ei (or, equivalently, 

is a switching of E 2 ). In words, it suffices to show that the probability that Ei is switched 
to E 2 is the same as the probability that E 2 is switched to Ei. To this end, we first construct a 
bijection (j) : Ei ^ E 2 between the edges of the two graphs, and then show that the conditioned 
probability measure is invariant under this bijection. The bijection is deterministic. 

Define E/\ := E 1 AE 2 and Ef-) ■= EiriE 2 , where A denotes the symmetric difference. Dehne 
IT := {// : e^(£'i) G £^a}- The interpretation of IT is the index set of neighbours of 1 in Ei that 
were switched in going from Ei to E 2 . This is the same set as the set from (I6.9p . Note that 
IT is now deterministic: for Ei and E 2 that are switchings of each other, the set IT is uniquely 
determined. Since E 2 is a switching of Ei, and by the constraints in the indicator functions 
and in the definition of Ts^s, we find that .EA|{ 2 ,...,iV} consists of |IT| disconnected subgraphs, 
each of which has 5 vertices and 4 edges. Each such subgraph is adjacent to a unique edge 
e^{Ei) where /r G IT. For /i G IT, we denote by the set of vertices consisting of 1 and the 
vertices of the subgraph that is adjacent to e^(£'i). By construction, E 2 \b^ is a switching of 
Ei\b^ (and vice versa); both are 1-regular graphs on six vertices. (The interpretation of is 
that the switching Ts,s that maps Ei to E 2 satisfies = Ei\b^-) This construction is illustrated 
in Figure EH 

Now we define the bijection </> : Ei —)• E 2 . For each /r G IT, we choose to be a bijection 
from Ei\b^ to E 2 \b^ such that 4‘{e^{Ei)) is incident to 1. (For each such /x there are two 
possible choices for this bijection; this choice is immaterial.) Without loss of generality, we can 
choose the enumeration e^{E 2 ) of E 2 such that e^(E 2 ) = (/>(e^(Ei)). This defines a bijection 
(j) '■ El Cl Ea —)• E 2 n Ea- We extend it to a bijection cj) : Ei ^ E 2 hy setting 4){e) ■= e for 
e G Er,. 

With these preparations, we now show (I6.12p . Given Ei and E 2 that are switchings of each 
other, we have constructed a set IT C [l,d] and subsets B^j^ for /x G IT, such that Ei\b^ and 
E 2 \b^ are 1-regular graphs obtained from a unique switching from each other. Since the are 
independent and since for each /x G IT the random variable is uniform on [1, 8], we find that 
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the left-hand side of (I6.12p is equal to 


8-1^1 p(5^ = ElI b, for^eff^, {/ie[l,dl:/^(Ei,S)J^(S) = l}=fy), (6.13) 

where S is uniform over 5i(Ei) x • • • x Sd{Ei). 

By an identical argument, we find that the right-hand side of (I6.12p is equal to 

8 -|H/| ^ ^^1^^ ioT GW, {fie [1, dj : 4(E2, S) J^(S) = l} = , (6.14) 

where S is uniform over Si{E 2 ) x • • • x Sd{E 2 ). Note that, by construction, W and for fj. G W 
are the same in both (|6.13p and (I6.14p . 

What remains is to show that ()6.13p and (16.141) are equal. In order to prove this, we 
abbreviate (j){S) := ... ,4>{Sd))- Then the dehnitions of /^, J^, and (p imply that for 

all S G Si{Ei) X • • • X 5rf(Ei) we have 

4(Ei,S) = 4((/)(Ei),(/)(S)), J^{S) = Jf,{cP{S)), (6.15) 

and hence (|6.13l) is equal to 

8-l'^lp(,/.(5^) = cPiE,)\B, for // G W, {/i G [1,^1 : /^(./.(Ei), ,/.(S)) J^(0(S)) = l} = w) . 

(6.16) 

Since (p{Ei) = E 2 and since 0 is a bijection from 5i(Ei) x ■ ■ ■ Sd{Ei) to Si{4>{Ei)) x ••• x 
Sd{4>{Ei)), and therefore (p{S) is uniform over 5 i(E 2 ) x • • • x Sd{E 2 ), we conclude that (I6.16P is 
equal to (I6.14|) . This concludes the proof. □ 

6.2. Estimate of exceptional configurations. In preparation of the proof of Proposition 13.71 for 
the uniform model, we now define the probability space H from Definition 13.51 The space 0 is 
the set of simple d-regular graphs (identified with their sets of edges E), and 

Ufj, ■= {sets of three distinct edges of the complete graph on N vertices} x |1,8] . 

Hence, the probability space Q = Q x Ui x ■ ■ ■ x Ud from Definition 13.51 consists of elements 

{e,ui,...,Ud) = (E, (S’i,si),...,(5d,Srf)) , 

where we denote elements of 0 by 0 = E and elements of by = (S"^, s^). Next, we define 
the (non-uniform) probability measure on the set D. To that end, we hrst endow the set of simple 
d-regular graphs 0 with the uniform probability measure. For each E G 0, we hx an arbitrary 
enumeration ei(E),... , erf(E) G E of the edges of E incident to 1. Then, conditioned on E G 0, 
we take tti,..., to be independent, with uniformly distributed over S^{E) defined in (j 6 . 6 p . 
and Sfj^ uniformly distributed over |1, 8 ]. (In other words, the probability measure is uniform on 
{6,ui,... ,Ud) G D such that each 5^ contains e^(E) and no other edge incident to 1.) Having 
dehned the probability space D, we augment it to D according to Dehnition 13.61 

From now on, we always condition on E and write = e^(E). We denote by p^, the two 
edges in not incident to 1 (ordered in an arbitrary fashion), so that 5^ = The 

next lemma provides some general properties of the random sets 5^. 

Lemma 6.2. The following holds for any fixed p G |l,(i]. 


35 




















(i) There are at most five v fJ- such that 

[5^] n [5j.] {1} and. [5i.] n [S^,] = {1} for all nfi . (6.17) 

(a) For any symmetric function F we have 

PiP^(l) + o(^^y\F\U- (6.18) 

Similarly, for any function F we have 

^G,J.,q^,F{pfJ,) = ]Eg^,p^F(g'^) = F(j>) + 0^—'^\\F\\oo ■ (6.19) 

Proof. We begin with (i). Let B be the set of u p satisfying (I6.17p . By definition, for n £ B, 
[S'jy] n [S'k] = {1} for all K Thus, each p G [S"^] \ {1} can be contained in at most one Si, 

with u £ B. The claim follows since [S^] \ {1} has at most five elements. 

Next, we prove (ii). Conditioned on the two edges p^^q^j, are by definition of chosen 
to be distinct and uniformly distributed on the Nd/2 — d edges not incident to 1. Let dl = {e £ 
-E : 1 € e} be the set of edges in E incident to 1. Then 

^gi,F{p^,q^) = Y. 

2 ) p,gGE\dl:p^q 

This shows (16.1811 : the proof of (I6.19jl is analogous. □ 

Next, we derive some basic estimates on the indicator functions and and the random set 
W. Ideally, we would like to ensure that with high probability IpJfi = 1. While this event does 
hold with high probability conditioned on Fq , it does not hold with high probability conditioned 
on Gfi. In fact, conditioned on it may happen that = 0 almost surely. This happens if 
there exists a, u p such that [5j^] n 7 ^ {1}. The latter event is clearly independent of To 
remedy this issue, we introduce the ^^-measurable indicator function 

h^ := l(e^nE^ = {l}), := IJ [S,] 

which indicates whether such a bad event takes place. Then, instead of showing that conditioned 
on we have IpJp = 1 with high probability, we show that conditioned on we have IpJp = h^ 
with high probability. This estimate will in fact be enough for our purposes, by a simple analysis 
of the two cases h^ = 1 and h^ = 0. In the former case, we are in the generic regime that allows 
us to perform the switching of with high probability, and in the latter case the switching of 
is trivial no matter the value of S^. 

For the statement of the following lemma, we recall the random set W defined in (j6.9p . and 
that is obtained from W by replacing the argument with We also recall that X/NY 
denotes the symmetric difference of the sets X and Y. 
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Lemma 6.3. The following holds for any fixed // G |1, dj. 

(i) Conditioned on Q^, with probability 1 — O(^), we have 

l(/iGTy) = = hf,. (6.21) 

Conditioned on Fq, with probability 1 — we have // G VL (i.e. h^ = l). 

(a) Almost surely, we have \W/\W^ \ {/i}| ^ 10. 

(Hi) Conditioned on Q^, with probability 1 — O(^), we have WAW^ = 0. 

Proof. We first show the claim of (i) concerning conditioning on Q ^. First, the definition of 
immediately implies that = 0 if = 0. Therefore and since /i^ is ^^-measurable, it 
suffices to show that, conditioned on such that h^ = 1, we have = 1 with probability 
1 — O(^). Hence, for the following argument we condition on and suppose that = 1. We 
estimate 

% iWn = 0 ) < % (4 = 0 ) + % (Ju = 0 ) ■ ( 6 - 22 ) 

The first term on the right-hand side of (I6.22p is equal to 

is not 1 -regular) 

= IPCm is not 1-regular) (6.23) 

+ 1-regular, is not 1-regular) (6.24) 

We first estimate (|6.23l) . Since p^ is uniformly distributed under the constraint p^ ^ 91, we 
find that (j6.23p is bounded by 0{d/{dN)) = 0{1/N). Similarly, given p^ such that Lile^^up^ 
is 1-regular, is uniformly distributed under the constraint ^ dl L) {p^}- Moreover, if 
H|e^up^ug^ is not 1 -regular then a vertex of must coincide with or be a neighbour of a vertex 
in Up^. From this we deduce that (|6.24l) is bounded by 0{d/{dN)) = 0{1/N). We have 
therefore estimated the first term on the right-hand side of (|6.22l) by 0{1/N). 

To estimate the second term on the right-hand side of (j6.22p . since 

{J^(S) = 0} = {[5^] n / {1}} C {e^ n / {1}} U {{p^ U q^) n / 0} (6.25) 

and since {e^ n = {1}} holds by assumption, it suffices to estimate lPg^((P/i U q^f) n 0). 
Clearly, \ {1} has at most 5{d — 1) vertices. This implies that p H 7 ^ 0 for at most 
0{d?) edges p. Taking F{p,q) := l(pn / 0) -|- l{q PI / 0) in (I6.18p . we therefore get 
(p^U( 7 ^)niL^ / 0 with probability 0{d/N). This proves that conditioned on with probability 
1 — 0(;^), the event (I6.2ip holds. 

Next, we show that conditioned on Fo, with probability 1 — 0{-^), we also have h^ = 1. By 
a union bound and since e^ Ci e^ = {1} if / p, 

^ J \ 

IPjFo(V = 0) = n / {1}) ^ ^PjFp(e^ n (pi. U O',,) / 0) = (6-26) 

as claimed. This completes the proof of (i). 

Next, we show (ii). We write W = and = S)). We then need to show |VFAW\{p}| ^ 
10. For this, we make the following observations. 
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1. If [5iy] n [Sk] / {1} for some v ^ n, with both distinct from then v regardless of 
Ufj,] therefore then also v and v ^ W /\W. 

Hence, all v G [W/\W) \ {/r} satisfy [5i,] H [S,^] = {1} for all n ^ {fi, i'}. 

2. Under this last condition, i.e. z/ / // satisfies [S'!,] n [5^] = {1} for all k ^ by 

definition G VU if and only if [5^] n [5^] = {1} and is 1-regular, and z^ G lU if and 

only if [S'!,] n [S^] = {1} and E\[s^] is 1-regular. 

Hence, i' G (lUAIU) \ {^} requires [5i.] fl [S^] ^ { 1 } or [^j,] n [Sfj] / {!}. 

We conclude that all v G {WAW)\{n} obey (I6.17P or (I6.17P with 5^ replaced by 5^. Therefore, 
Lemma [Ql il implies |(IUAIU) \ {^}| ^ 5 -|- 5, as claimed. 

Finally, we establish (iii). For this, observe that lu is independent of and that Jy is 

independent of Pij,,q^ on the event {(p^ U q^) PI = 0}. In the proof of (i), we have already 
shown that the latter event has probability at least l — conditioned on Qfj_. This concludes 

the proof. □ 

6.3. Proof of Proposition 13.71 for the uniform model. With the preparations provided in Sec¬ 
tions 16.1116.21 we now verify the claims of Proposition 13.71 for the uniform model. 

Proof of Proposition [3771 uniform model. The parametrization obeying Definitions [33L 
13.61 was defined at the beginning of Section [6121 The random variables ai,..., a^, ai,..., a^, and 
A are defined as follows. By definition, A is the graph with edge set Ts,siE), i.e. 

A := M{Ts4E)). (6.27) 


Moreover, is by definition the unique vertex incident to 1 in the subgraph 


rs,s(^)|[5,] if/iGiu 
{e^E)} if/i^W, 


where we recall that W was defined in (16.9p . The definition of is illustrated in Figure [631 
Then by Lemma l 6 .ll A is the adjacency matrix of the uniform model. To see that ai ,... ,ad 
are an enumeration of the neighbours of 1, it suffices to show that ^ for p ^ v. This 
follows from the following simple observations: if p, u ^ VF then if p, u G VF then by 

dehnition of W we have [5^] n [5iy] = {1}; if p G VF and v then by dehnition of W we have 
&v n [Sfj] = {1}. This proves (i). 

What remains is the definition of and the proof of (ii) and (iii). To define o^, we denote 
by p^ and the two edges in not incident to 1, ordered in an arbitrary fashion. (Note that, 
by definition of 5^ from (j 6 . 6 p . e^, p^, and q^ are distinct but not necessarily disjoint.) We label 
the vertices of = {p^,p^} and = {qj^, q'^} in an arbitrary fashion, and take the pair (a^, b^) 
to be uniformly distributed (parametrized by G |1,8]) in the set 


{(pI, qI), (pI, qI), (pI^qI), (qI^pI), (qI^pI), (qI^pI), (qI^pI)} ■ (6-28) 


More precisely, we parametrize (a^, b^) = ( 0 ^( 5 "^, s^), bfj,{Sfj,, s^)), with G |1, 8 ], in such a way 
that, in the nontrivial case = 1 , the switching T 5 ^^s^(£') from (16.2h is given as in dSS]) by 


M{Ts, 


,,,s^{E)) — Tr^l,a^a b^b 


(Mm, 


(6.29) 
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^2 


Figure 6.4. Solid lines depict edges oi E\ Ts,s{E), dashed lines edges of Ts,s{E) \ E, and dotted lines 
edges of F n Ts,s{E). For the graph E, the two encircled groups contain the three edges of Si and the 
three edges of S 2 respectively. In the switching E >->• Ts,s{E), the solid edges are switched to the 
dashed conhguration; hence, 1 G W and ai is the unique neighbour of 1 in [S'!] after the switching. On 
the other hand, S 2 is not switchable (so that 2 ^ W) since |[5'2]| < 6, and 02 is given by the original 
neighbour of 1 in 62 . 


where r is defined in ()3.3I) . and the vertices are defined by the conditions = {l,r^} 

and (/^} = {{o^, a^}, Note that if are disjoint, is uniformly distributed 

on U and 6^ uniformly distributed on p^ or whichever does not belong to. 

We shall show (ii) and (iii) with the high-probability events given by those on which the 
conclusions of Lemma 16.31 hold. More precisely, the high-probability event in (iii)(l) is given by 

= K} n = 0} (6.30) 

and the high-probability event in (ii)(2) and (iii)(2) is given by 
S;. := { 4 = 1} n = 0} = {/r G W} n {WAW^ = 0} = H^n{V = l}. (6.31) 

By Lemma [6]3] (i) and (iii), recalling that Fq C G^i and ^ by (II.5p . 

We now show (ii). Prom the definition of (a^, 6^), we find that is chosen uniformly among 
the four (not necessarily distinct) vertices pji,pfi, qji, gfi- Therefore we get, for any function / on 
[1, N] that 

1 ^ /I 

i=l ^ 

where we used (|6.18p with F{pfj,, q^) := |(/(p^) + fipfi) + fiqji) + f{qfi)) and the fact that each 
vertex is contained in exactly d edges, so that I]{ei,e2}e£; + /(^^)) = Tr /(*)• 

This shows (ii)(l). To verify (ii)(2), we use that on the event we have p G W. Moreover, from 
(I6.29P we find that on we have {l,o^} G Ts^^sij,{E), and consequently, using the definition of 
W and (I6.10p . {l,a^} G Ts^s{E). Since G [S^] the definition of immediately implies that 
Q/i = 0 (^ on B^. Together with (I6.32p . this concludes the proof of (ii)(2). 

To show (iii), we often drop the sub- and superscripts p and abbreviate a = A = 
and so on. For the first claim of (iii)(l), it suffices to show that ^ — A is always a sum of at 
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most 72 terms 'iz/S.xy By (16.111) 


A-A = W){M{Ts^^s,{E)) - M{E)) - e W){M{Ts^ -^{E)) - M{E)) 

+ ±{M{Ts^,s^m-M{E)), (6.33) 

u£WAW\{ij.} 


(where the sign ± in the last sum is + if £ hh and — if £ hh). In Lemma [6.31 (ii), we proved 
that |ILAiy\{/r}| ^ 10. Therefore, since each term M {Ts^s (E)) — M{E) is the sum of six terms 
±Aa;j^ by (|3.3p . we find that ^4 — ^4 is the sum of at most 12 x 6 = 72 terms EA^y, as desired. 
This proves the first claim of (iii)(l). 

Next, we verify the second claim of (iii)(l) and (iii)(2). By (|6.32p we may assume that the 
event Ey holds. In particular, since IT AIT = 0, (16.lip implies 

71 - i = h[M{Ts^,s,m - M{Ts^ -^iE))) . (6.34) 

In the case h = 0, the right-hand side vanishes and the second claim of (iii)(l) is trivial. On the 
other hand, if /i = 1, by (j6.29p . 


A-A = M{Ts^,s^{E)) - M{Ts^^,^{E)) = Tri,aa,i,k{M{E)) 


— T 


rl,aa,66 


{M{E)). (6.35) 


As in the proof of (ii)(l), we find that conditioned on each of the random variables a, b,a,b 
is approximately uniform, and using (I6.19jl instead of (I6.18jl . that b,b are each approximately 
uniform conditioned on Qy and a, a. The same holds with a,b,a,b replaced by d,b,a,b. Thus, 
under the probability distribution conditioned on Gy and the event {h = 1}, we have 


A — A — Ala ~ Ala + A , 


(6.36) 


where X is of the form (|3.8I) and has the property (a) or (b) from Remark 13.31 (with all prob¬ 
abilities given by conditioning on Qy). Since for each term EAxy in X at least one of x and 
y is approximately uniform, and since a and d are approximately uniform, we conclude that, 
conditioned on Qy and the event Ey D {h = 1}, A — A is given by a sum of 10 terms EAxy 
such that at least one of x and y is approximately uniform. Together with the trivial identity 
A — A = 0if/i = 0 and (I6.32p . this shows the second claim of (iii)(l). 

Since T,y = Ey f] {h = 1}, the identity (I6.36P also holds on the event Recalling (I6.32p 
and the form of X from (13.811 . we obtain (iii)(2). This concludes the proof. □ 


7. Permutation model 

In this section we prove Proposition 13.71 for the permutation model. First, we consider the 2- 
regular random graph defined by a single uniform permutation. As before, the symmetric group 
of order N is denoted by S'at, and for any permutation a G Sn, the associated symmetrized 
permutation matrix is denoted by 


N 

P{a) = P{a-^) := 

i=l 


(7.1) 
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b+ f tt 1(6+) 

^ n/ 

/ \ 

a+ ^-• 7r“i(a+) 

/ ^ s 

2^ ->1 

a_ •—\ ^ 7r(a_) 


b- •---7r(6_) 


6+ ^ TT 1(6+) 

2 :: ji 1 = a+ 

^ / 
y 

/ \ 

a_ •-\7r(a_) 

\/ ^ 

6_ •-- '^{b-) 


Figure 7.1. The solid arrows depict the permutation tt, and the dashed arrows the permutation 
'7a+,a_,f)+,f)_ (tt) defined in (17.3L The left diagram depicts the generic case in which all of the elements 
of {l,2,a+,a-, 6+,6_,7r i(a+),7r i(6+),7r(a_),7r(6_)} are distinct. The right diagram depicts the case 
where a+ = 1 but the remaining vertices are distinct; this leads to a loop at 1. 


Denote by 'jij = jji € Sn the transposition that exchanges i and j. For the remainder of 
this subsection, we identify Sn -2 as the subset of Sn of permutations that exchange 1 and 2, 

Sn -2 = {ttgSn-- 7r(l) = 2,7r(2) = 1} . (7.2) 

For TT G Sj\f -2 and a+,a_,6+,6_ G |l,iV|, we define Ta^a-b+b-iT^) £ Sn by 

Fa+a_b+6_(vr) := l2b+l2a+T^l2a-l2b. ■ (7.3) 

As illustrated in Figure [7T1 in the case that {1, 2, a+, a_, 6+, 6_, 7r“i(a+), 7r“i(6+), 7r(a_), 7r(6_)} 
are distinct, the action of T on vr amounts to two double switchings, as depicted in Figure ItTI 

Lemma 7.1. //(tt, a+,a_, 6+, 6_) is uniform on Sn-2 x x |2,Aip, then Ta^a-b+b-iT^) is 

uniform on Sn- Moreover, 

(Fa+a_6+fe_(vr))(l) = a+, iTa+a_b+b_{T^))~^{i) = , (7.4) 

provided that |{1, 2, o+, a_, 6+, 6_}| =6. 

Proof of Lemma [7.11 First, we show that the map a : Sn -2 x |1, A^l x |2, A^] — >• Sn defined 
by 

(7r,a+,a_) i —a = d{7T,a+,a-) := 72a+7r72a_ (7.5) 

is a bijection. This follows from the following explicit inverse map d i—)• (7r,a+,a_): 

(i) if d(l) / 1 then a+ := d(l), a_ := d“^(l), and vr := 72 a+o' 72 a_; 

(ii) if 5'(1) = 1 then a+ := 1, a_ := d“^(2), and vr := d'j 2 a-- 
Therefore, for all fixed 6 +, 6 - G |1, A^], also the map 

S'Ar _2 X |1, A^] X |2, A^] —^ Sn, (7r,a+,a_) i—)■ 726+72a+7r72a_72fe_ (7.6) 
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is a bijection. In particular, under (17.61) . the uniform distribution on Sn -2 x x [2, A^] is 

pushed forward to the uniform distribution on Sj\f. Clearly, the distribution remains uniform 
after averaging over the independent random variables 6+,5_ G |2,A^]. Finally, (17.41) is easily 
verified. This completes the proof. □ 

The probability space for the permutation model (see Definition 13.5p is realized as (13.171) 
endowed with the uniform probability measure, where 0 := and 


:= [l,iV| X [2, AT] 

(/zG [l,(i/2|). 

(7.7) 

:= [2,iV| X [2, AT] 

(m £ 1^/2 + 1, rf]) ■ 

(7.8) 

Elements of 0 and are written as 9 = (vri 

,..., TTd/ 2 ) e 0 and = {a^, b^) 

G U^. Eor 

/i G [1, (f/2] we define the random variable 





(7.9) 


By Lemma l7.ll ui,..., (t ^/2 i.i.d. uniform permutations in Sjsf- The adjacency matrix of the 
permutation model is given by 

d/2 

A ■.= ( 7 - 10 ) 

ti=i 

It is convenient to augment the sequence (cr^) to be indexed by |1, dj by defining := for 
/r G ld/2 + 1,/i]. Hence P{(yfi) = P{(Td-^) for all /r G |l,(i]. Also := where /r G 

is an enumeration of the neighbours of 1 in A. 


Proof of Proposition 13.71 permutation model. We use the parametrization of the prob¬ 
ability space defined above, which satishes the conditions of Definition 13.51 and augment it 
according to Definition 13.61 Then claim (i) follows immediately from Lemma l7.ll 

To show (ii), we hrst recall that is uniform on |1, A^] if /r G |1, (i/2] and that is uniform 
on |2,A^| if /i G |(i/2-|- 1, (i|. Either way, conditioned on is approximately uniform on 

|1,A^|. Moreover, for fj, G |l,(i], conditioned on Pq, with probability 1 — 0{^), the event 
{\{l,2,afj,,b^,ad-f_t,bd-fi}\ = 6} holds. By Lemma ITTI on this event we have, for fi G |l,(i/2|, 
= <7^{l) = and ad-fj, = This concludes the proof of (ii). 

What remains is the proof of (iii). We fix /i G |1, (i|, and drop the index ^ from the notation 
and write A = a = and (t_ = CTd-^i-, and so forth. Then set 


Q ■= 72&_72a_7r 

so that a = ^ 2 h'~l 2 aQ~^■ Given set Ag{a,b) ■= P{'^ 2 bl 2 aQ~^)- Then, for all a,b,d,b, 

Ag(d, b) — Ag{a, b) = Agd,b — Aga,b + Aga, 726722 “ ~ ^£'“,726720“ + Ag 2 ,'Y 2 i,a ~ ^^ 2 , 725 “ 
Ag{a, b) — Ag{a, b) = 2 ~ At-;,^2 + + ^r2,b ~ AT- 2 ,fe , 

where r := Q^ 2 a- Indeed, to verify (17.121) . note that Ag{d,b) — A^{a,b) is given by 


N 


n 

i=l 


N 


7 ,7267225 Td ^d72672a5 Md 


= 


i=l 


5(d.726722(d ^5(d.72672a(d ) ’ 


(7.11) 


(7.12) 

(7.13) 


(7.14) 
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and that the differences under the last sum vanish unless i G {2, a, o}, and are given by (17.121) . 
Similarly, for the difference Ag{a, b) — Ag{a, b), we obtain 

N N 

(i) ~ ‘^*,72i,72Qe“bd) ~ “ ^T(*).72b(*)) ’ 

i=l i=l 

and observe from this representation that only i G {2, b, b} yields a nonzero contribution, and 
that the corresponding terms are given by (j7.12j) - (l7.1.Sp . 

To prove (iii), observe that only the term P(cr^) = P{ad-fi) in (I7.10p contributes to A—A and 
that therefore A —A = ^^(a, b) — Ag{a, b). Prom (I7.12jl - (j7.13ll . it is straightforward to verify (iii). 
The first statement of (iii)(l) holds since (j7.12p - (l7.13p contains at most 12 terms P/S.xy. For the 
second statement of (iii)(l), note that, since g is independent of a,a,b,b, and r independent of 
b,b, for each of these terms Axy in (I7.12l) - (j7.13p . at least one of x,y is approximately uniform 
since a, a, b, b are approximately uniform. For (iii)(2), observe that, conditioned on Tq, £»^(2) = 1 
and 2, a, d,5, 6 are distinct with probability 1 — 0{^). On this event, (iii)(2) can be verified 
directly from (I7.12p - (j7.13p . For example, A^ 2 , 726 a ~ ^s‘^,'y 2 ba = ^la in ()7.12p . We skip 

the details for the other terms. □ 


8. Isotropic local law and probabilistic local quantum unique ergodicity 

In this section we state and prove the isotropic local semicircle law for A, which controls the 
difference G — ml in the sense of generalized matrix entries (a,Gb) — m(a,b), instead of the 
standard matrix entries from Theorem 11.11 obtained by taking a and b to lie in the standard 
coordinate directions. The arguments used in this section rely crucially on the exchangeability 
of the random regular graph. This is different from the remainder of the paper, in which we 
did not exploit exchangeability in an essential way. An isotropic local law was first proved for 
Wigner matrices in [32] and subsequently extended to generalized Wigner matrices and sample 
covariance matrices with uncorrelated population in [5] . Recently, such control was also obtained 
for sample covariance matrices with general population [33], in which case G is approximated 
by an anisotropic matrix that is not a multiple of the identity. 

As applications of the isotropic local law, we establish the isotropic delocalization of eigen¬ 
vectors ('Corollary 18.4p and a local quantum unique ergodicity result fCorollary 18.51) . In the 
following we call an ^^-normalized vector a unit vector, and write a T e if ^ • ajCj = 0. 

Theorem 8.1 (Isotropic local law for random regular graphs). Under the assump¬ 
tions of Theorem M.il for any deterministic unit vectors a, b T e, any S> 1, and any z G C+ 
satisfying rj (,‘^/N, we have 

(a,G( 2 :)b)-m(z)(a,b) = 0(Tj(^^>(z)) + ^C^^>( 2 :)) (8.1) 

with probability at least 1 — . 

Note that the isotropic law in the subspace spanned by e is trivial since G{z)e = —z~^e. 
Theorem 18.11 follows immediately from Theorem 11.11 and the following general result for ex¬ 
changeable random matrices. Recall that a random vector (I^)(^^ G is called exchangeable 
if for any permutation a £ Sn have 

(y,), ^ (y.(,)),. (8.2) 
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Similarly, a random matrix G £^NxN jg called exchangeable if for any a G Sn we have 

0^ij)i,j = ■ (8-3) 

In particular, the (normalized and centred) adjacency matrices of any of the models of random 
d-regular graphs introduced in Section [1.21 are exchangeable. 

Theorem 8.2 (General isotropic local law). Let G be the Green’s function (11.81) of an 
exchangeable random matrix H at some z G C+ . Then, for any deterministic To, > 0, 
m G C, unit vectors a, b T e, and C 1, we have 

(a,Gb) - m(a,b) = 0(Trf + C^To) (8.4) 

with probability at least P(maxj|Gjj — m| ^ T^, maxj^jjGjjl ^ Tq) — 

The proof of Theorem 18.21 follows from the following moment bounds for exchangeable ran¬ 
dom matrices. The estimate (18.5|) was previously established for p = 2,4 in HIE]. 

Proposition 8.3. Let ai, ... ,aN G C 6e deterministic with = 0 and ^ 1- 

(i) Let be a exchangeable random vector. Then for allp'^1 we have 

N 

i=l 



(a) Let (Tjj)^-^^ be a exchangeable random matrix. Then for all p ^ 1 we have 


N 

diO-jYij 

*J=1 





( 8 . 6 ) 


The proof of Proposition 18.31 is given in Appendix [Bj 

Proof of Theorem 18.21 By polarization and homogeneity, it suffices to consider the case 
where a = b is a unit vector perpendicular to e. Define Yij := (/)l(i 7^ j)Gij with the indi¬ 
cator function cf := 1 (maxj^jjGjjj ^ Tq). Then {Yij) is an exchangeable random matrix. By 
Proposition 18.31 (ii) with p = CVlog C C S> 1 , we get using using Markov’s inequality 


Ka,ya)| >CC%c 



for some constant G. Since for any unit vector a we have (/>|(a, [G — m)a)| ^ max* \Gii — 

|(a, Ta)!, the claim now follows by a union bound. □ 

The isotropic local law implies the isotropic delocalization of the eigenvectors of A, which 
also follows from Corollary 11.21 and Proposition 18.31 lib similarly to the proof of Theorem 18.21 

Corollary 8.4 (Isotropic eigenvector delocalization). Under the assumptions of The- 
orem{LJ[ for any unit eigenvector v of A or H, any deterministic unit vector a T e, and any 
1, we have (a, v) = 0{f,Cf^/y/N) with probability at least 1 — 
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Finally, we note that Corollary 11.21 and Proposition 18.31 (i) also imply the following local 
quantum unique ergodicity result, similarly to the proof of Theorem 18.21 

Corollary 8.5 (Probabilistic local quantum unique ergodicity). Let a : |1, A^] —)■ M 
be a deterministic function satisfying = 0- Under the assumptions of Theorem \l.ll for 

any unit eigenvector v = {vi)fLi of A or H and for any 1, we have 

£«<»? = (8.7) 

with probability at least 1 — 

Corollary [83] states that, on deterministic sets of at least vertices, all eigenvectors of the 

random graph A are completely flat with high probability. In other words, with high probability, 
the random probability measure i vf is close (when tested against deterministic test functions) 
to the uniform probability measure i i—)• 1/N on [1,A^]. For instance, let I C |1,A^] be a 
deterministic subset of vertices. Setting a* ■■= l{i G I) — \I\/N in Corollary 18.51 we obtain 





JV ) 


( 8 . 8 ) 


with probability at least 1 — — e“^. The main term on the right-hand side of (18.8p is 

much larger than the error term provided that |/| (^^C)^- Note that we can obtain (18.711 and 

(18.81) asymptotically almost surely with (CC)^ = (logA^)^ by choosing ^log^ = C(logA^)^ and 
= C~^ log^ for some large enough constant C > 0, so that (18.8p is a nontrivial statement for 
|/| ^ (logiV)^ 

The celebrated quantum chaos conjecture states that the eigenvalue statistics of the quanti¬ 
zation of a chaotic classical system are governed by random matrix theory [61 140p4H l46] . Random 
regular graphs are considered a good paradigm for probing quantum chaos; see [l2] for a review. 
For generalized Wigner matrices, a probabilistic version of QUE, as well as the Gaussian distri¬ 
bution of eigenvector components, was proved in m- The first result of this kind, for a smaller 
class of Wigner matrices, was obtained in |31ll43j . Moreover, using the local law proved in this 
paper, these results can be extended to the d-regular graph as well [TO] . 


Remark 8.6 (Erdos-Renyi graphs). We conclude this section by remarking that all of the results 
from this section - Theorems 18.11 and 18.21 and Corollaries 18.41 and 18.51 - have analogues for the 
adjacency matrix of the Erdos-Renyi graph, whose proofs follow in exactly the same way using 
Proposition 18.31 and [20l Theorem 2.9]. We leave the details to the interested reader. 


A. Improved bound near the edges 


In this appendix, we sketch the changes required to improve (I1.12p such that is replaced by 
(|1.20p on the right-hand sides, namely by 





where T 



1 

7B- 


(A.l) 
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First, analogously to F and F^, define 

T := maximGii, := ||T||£,oo(g^) . (A.2) 

Then it is easy to see that Lemma 13.91 (ii) can be improved to 

T. 1 F2 

— H-^- — . 

Nr] Nr] /d 

Assuming that F^ = 0(1) and that = 0(5), we have = 

5 1 1 

Nr]~^ Nr]^ y/D ' 

Next, Proposition 12.21 can be improved as follows. 

Proposition A.l (Improved version of Proposition 12.21b Suppose that ^ > 0, C > 0, and 
that O S> Let 6 ^ N~^ he deterministic. If for z we have 




replace F^<I>^ by <I>^ where 

(A.3) 

0 (<I> 5 ) where 

(A.4) 


r*{z) = 0(1), T ^ 5 
with probability at least 1 — e“^, then 

max.\Gii{z) - m{z)\ = 0{F^{C^5{z))), max|Gjj(^;)| = 0(^4>5(^;)), (A.5) 

* 

with probability at least 1 — e“(Oog?)AC+0(iogA^)_ 

Sketch of proof. We first verify that in all estimates of Sections HHSl the parameter <I> arises 
from just two possible sources: from or Lemma 13.91 (ii). 

In particular, by the improved version of Lemma [T9] discussed around (IA.3p . Lemma [4. 6 1 can 
be improved so that F® is replaced by F^<I>^. This implies an improved version of Proposi¬ 
tion [TT] in which the assumption that F = 0(1) holds with probability at least 1 —e“^ is replaced 
by the assumption F = 0(1) and T = 0(5) with the same probability, and <I> is replaced by 
in the conclusion. The proof of the improved version of Proposition 14.II is then analogous to the 
proof of Proposition 14.11 given in Section |3| In particular, note that using 5 ^ N~^, it is easy 
to verify that there are Yjj satisfying all required conditions and F^^$^ ^ 

Similarly, given the improved versions of Lemma 13.91 and Proposition 14.11 the proof of the 
improved version of Proposition 12.21 is identical to that given in Section (5] □ 


Finally, given Proposition lA.ll the proof of Theorem 11.11 involves a slightly more involved 
induction than that given in Section [2l in which we propagate both estimates 


F ^ 2, 


T < 


^Imm -|- 



(A.6) 


simultaneously, for some sufficiently large constant Q. The hypothesis ()A.6p is trivial for p ^ 1. 
By an explicit spectral decomposition (or by an argument analogous to the proof of Lemma l2.1l) . 
it is easy to verify that p Im Gu is increasing in r], and hence (IA.6P for some z = E + ir] implies 


T{z’) ^ 5(z') 


4(5^ Imm(z') -|- \ ^'I'(z') -|- 


Nr]' 


2/Y 


(A.7) 
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for z' = E -\- \r]' with rj' = t]/2. Hence, we may apply Proposition lA. 11 It is easy to verify that 


f-I-i < (A.8) 

Therefore, using Fz{r) ^ y/r, we get from Proposition lA.II that 

T{z') ^ Imm{z') + 0{y/^s) ^ 0(Q) ^Imm+ (A-9) 

and this propagates the induction hypothesis (IA.6p since 0{Q) ^ for a sufficiently large Q. 


B. Moment ho un ds for exchangeable random matrices: proof of Proposition IST^ 

In this appendix we prove Proposition 18.31 To avoid extraneous notational complications arising 
from complex conjugates, we suppose that all quantities are real-valued. We abbreviate [n] := 
[1, n], and denote by the set of partitions of |n]. By Holder’s inequality, it suffices to consider 
the case p £ 2N. 

We begin with (i). Abbreviate X := For i £ {Nj^ define P(i) £ as the partition 

generated by the equivalence relation k ^ I and only if ik = ii- Then we get 

= E = E ^ l(P(i) = H) n , (B.l) 

ieIAf]pfc=i k=i neTp iel^vp k=i 

where A'(n) := lij, for any i satisfying P(i) = H. That Kill) is well defined, i.e. 

independent of the choice of i, follows from the exchangeability of (Yi). For future use we also 
note that by Holder’s inequality we have |A'(n)| ^ ||Fi||p. We use the notation vr £ H for the 
blocks of H. Next, we rewrite the sum over i £ {NJP as a sum over r = (?' 7 r) 7 ren G {xj^ to get 

* 

^ iF(n) H 4:' , (B.2) 

negJp reliv]n7ren 

where the star on top of the of the sum indicates summation over distinct indices r, i.e. 

* 

E = E If (b,3) 

r r ee£(n) 

where T(n) := {{vr, vr'} : vr, vr' £ H, vr / vr'} is the set of edges of the complete graph on the 
vertex set H, and := 1(7V = ^Tr')- 

We need to estimate the right-hand side of (jB.2p by exploiting the condition Oj = 0. 
Obtaining the bound of order (p^/logp)^ requires some care in handling the combinatorics. 
We shall multiply out the product in (IB.31) . which has to be done with moderation to avoid 
overexpanding, since the resulting sum is highly oscillatory. The naive expansion n ee£:(n)(l 
le) = X]£;c£:(n) We&Ei~^e) rough. Instead, we only expand a subset of the edges ^(H), 

and leave some edges e unexpanded, meaning that the associated factors (I — le) remain. 
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The partial expansion of the product nee£’(n)(^ “ is best formulated using edge-coloured 
graphs. We consider graphs on 11 whose edges are coloured black or white, and denote by B 
the set of black edges and by W the set of white edges. Thus, an edge-coloured graph is a pair 
{B, W) C T(n)^ satisfying i? n IT = 0. For any edge-coloured graph {B, W) we define 

Jb,w ■■= Hi-Qllil-Ie). (B.4) 

eSS e&W 

Hence, each black edge e £ B encodes the indicator function —If, and each white edge e G W 
the indicator function 1 — Ig. Note that for e € VF we have the trivial identity 

’^B,w = -^B,W\{e} + ’^BU{e},W\{e} ■ (B.5) 

We shall define a process on the set of edge-coloured graphs that operates on each white edge, 
either leaving it as it is or generating two new graphs using (IB.5p . one graph where this white 
edge is removed, and another graph where the white edge is replaced by a black one. To that 
end, we choose a total order on T(n) and denote by e— and e+ the immediate predecessor and 
successor of e. We denote by emin and Cmax the smallest and largest edges of T(n), and introduce 
the formal additional edge 0 to be the immediate predecessor of emin- 

For each e G T(n) we shall define a set of ^n(e) of edge-coloured graphs {B,W) such that 
W contains all edges greater than e. The sets Gn(^) are defined recursively as follows. First, we 
set 0n(O) := {(0,T(n))}. Thus, ^n(O) consists of the complete graph with all edges coloured 
white. Then for emin ^ e ^ Cmax the set Qn{e) is obtained from 0n(e—) by 

Gn{e) ■■= U U{B,W,e), (B.6) 

(-B,VK)een(e-) 

where lJ{B,W,e) is a set of one or two edge-coloured graphs obtained from {B,W), using one 
of the two formulas 


ll{B,W,e) = {{B,W\{e}),{BU{e},W\{e})}, (B.7) 

l({B,W,e) = {{B,W)}; (B.8) 

which choice among ()B.7p and (jB.Sp to make will be determined in (jB.lOp below. The choice 
()B.7p amounts to multiplying out 1 — Ig and (jB.Sp to not multiplying out 1 — Ig. Note that, by 
construction, we always have e C PF on the right-hand side of (IB.61) . Moreover, by (IB.Sp . no 
matter which choice we make between (IB.7p and (jB.Sp . we always have the identity 

Jb,w = ^ Jb,w 

{B,W)&gu{e-) {B,W)&gn{e) 

for all e G T(n), and hence by induction 

n (l-/e) = Jb,W- (B.9) 

eeS{n) (B,M/)een(emax) 

Note that always choosing ()B.7p leads to the identity nee£^(n)(^ ~ ^e) = 'YjB<zS{n)Y\.eeBi~^(i)-, 
which, as explained above, is too rough; conversely, always using (IB.81) leads to the trivial 
identity neef(n)(l “ ^e) = nee£(n)(l “ h)- 
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Figure B.l. The process from C/n(0) to t/n(emax)- Black and white vertices are depicted using black and 
white dots, respectively. Black and white edges are depicted using solid and dotted lines, respectively. 
Since the complete graph on 5 vertices has 10 edges, there are 10 steps, i.e. maps from a coloured graph 
to another. We start from the complete graph whose edges are all white. We highlight the white edge 
e indexing the subsequent step by drawing it using dashed lines. In each step using (jB.7|) . we choose 
to draw one of the two possible resulting graphs. The steps 6,8,9,10 use (IB.8I1 . and the other steps use 

(EJl). 


In order to define which choice of U in (IB.7ll - (jB.8h we make, we also colour the vertices 11 
black or white. A vertex is black if it is a block of size one and white if it is a block of size 
greater than one. This defines a splitting of the vertices 11 = Hi U 112 into black vertices Hi 
are white vertices 112, and also induces a splitting of the edges ^(n) = Ti(n) U Ti 2 (n) U ^2(11), 
where £iij\) is the set of edges connecting two vertices of Ilj for i = 1,2, and Ti 2 (n) is the set 
of edges connecting two vertices of different colours. We choose the total order on T(n) so that 
Ti(n) < Ti 2 (n) < ^ 2 ( 11 ). With this order, we define our choice oiU: 


use 



if e is incident to a black vertex that is not incident to a black edge, 
otherwise. 


(B.IO) 


See Figure [BT] for an illustration of the resulting process on coloured graphs. 

Let {B, W) G ^n(emax)- The following properties can be checked in a straightforward manner 
by induction: (a) W is uniquely determined by B (given the colouring of the vertices and the 
total order on T(n)); (b) i? is a forest (i.e. a disjoint union of trees); (c) a black vertex can only 
be incident to a white edge if it is also incident to a black edge; (d) two white vertices cannot 
be connected by a black edge; (e) a black and a white vertex can only be connected by a black 
edge if the black vertex is not incident to any other black edge. 

Now going back to (lB.2h . we find using (IB.91) that 


^ K{U) E '^sw(r) n 4:' (B.ll) 

neTp (B,VU)een(emax)re[Ar]n ttSH 

By property (c), if there is a black vertex that is not incident to a black edge, it is also not 

incident to a white edge, and Oj = 0 therefore implies X^re[Af]n'^B,vu(r) flTren ~ 0. 
Therefore the sum over [B, W) can be restricted to graphs in which every black vertex is incident 
to at least one black edge. For such graphs, ^ ^ \Jb,w\ ^ 1 imply the bound 
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|Sre[Ar]n <^B,VK(r) riTrenoil'l ^ 1- Using |K(n)| ^ ||yL||fJ we conclude 

EX^ ^ ||yi||P J]|gn(e„,ax)|. (B.12) 

ne<Pp 

It remains to estimate the sum on the right-hand side of ()B.12p . For hxed H, by (a) above it 
suffices to estimate the number of B satisfying the remaining conditions (b)-(e). From now on, 
all graph-theoretic notions always pertain to B, i.e. we discard all white edges. Let denote 
the set of black vertices not adjacent (by a black edge) to a white vertex: 

$ := {tt S Hi : vr is not adjacent to a vertex of 112 } . 


We shall estimate the number of graphs B on FI associated with any fixed ‘h. By (e), each 
vertex tt G Hi \ <I> has degree at most one. With (d), this gives the upper bound llli \ qj^ 

the possible choices of B in 11 \ <I>. Moreover, by (b), B is a forest on $. By Cayley’s formula 
^n -2 ^ number of trees on n vertices and the bound ^ (n/logn)"' on the 

number of partitions of a set, we find that there are at most (|^p/log|$|)l^l forests on <h. In 
summary, we conclude that the number of graphs B associated with 11 and C Hi is bounded 
by (|$|2/log|$|)l'^l(|ni| - |^>|)|n 2 l. 

Abbreviating k = |ni| and / = |d>|, we therefore obtain 



for some universal constant C > 0, where the factor (|) accounts for the choice of IIi, the factor 
{{p ~ k)/iog{p — k)y~^ for the choice of 112 , the factor (^) for the choice of and the factor 
{P/\ogiy{k — l)P~^ for the choice of B as explained above. Here in the last inequality we used 


p — k 
log(p — k) 


p—k 

{k-iy-^ ^ 





k 


for 0 ^ ^ /c ^ p. This concludes the proof of (i). 

Next, we prove (ii). By splitting Y into its diagonal and off-diagonal entries and using 
Minkowski’s inequality, it suffices to prove (j8.6p under the assumption = 0 for all i. Similarly 
to the proof of (i), we write 



2p p 


ieIA^]2pfc=i k=l 


2p 

Y, KiU) Y i(ni) = n)n«*.> 


neT2p ielAffp k=i 


where K{IV) := '^hk-ihk ^ satisfying P(i) = H (recall the definition of P(i) G ^ 2 p 

above (IB.ip i. As in the proof of (i), Ar(n) is well-defined by exchangeability of (Yij). By 
Holder’s inequality and exchangeability, we have the bound |iF(n)| ^ ||yi2||p- Now the proof 
of (i) following (IB. lb may be taken over verbatim, by replacing p with 2p. This concludes the 
proof of (ii). The proof of Proposition 18..SI is therefore complete. 


50 















Acknowledgements 


AK was partly supported by Swiss National Science Foundation grant 144662. HTY was partly 
supported by NSF grants DMS-1307444 and DMS-1606305. HTY and RB were partly supported 
by a Simons Investigator Award. The authors gratefully acknowledge the hospitality and support 
of the Institute for Advanced Study in Princeton, and the National Center for Theoretical 
Sciences and the National Taiwan University in Taipei, where part of this research was carried 
out. The authors’ stay at the IAS was supported by NSF grant DMS-1128155. 


References 

[1] N. Alon. Eigenvalues and expanders. Combinatorica, 6(2):83-96, 1986. Theory of computing 
(Singer Island, Fla., 1984). 

[2] N. Anantharaman and E.L. Masson. Quantum ergodicity on large regular graphs, 2013. 
Preprint, arXiv; 1304.4343. 

[3] R. Bauerschmidt, J. Huang, A. Knowles, and H.-T. Yau. Bulk eigenvalue statistics for 
random regular graphs. Ann. Probab., 2016+. To appear. 

[4] E. Benaych-Georges and A. Guionnet. Central limit theorem for eigenvectors of heavy tailed 
matrices. Electron. J. Probab., 19:no. 54, 27, 2014. 

[5] A. Bloemendal, L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Isotropic local laws for sample 
covariance and generalized Wigner matrices. Electron. J. Probab., 19:no. 33, 53, 2014. 

[6] O. Bohigas, M.J. Giannoni, and C. Schmit. Characterization of chaotic quantum spectra 
and universality of level fluctuation laws. Phys. Rev. Lett., 52:1-4, Jan 1984. 

[7] C. Bordenave. A new proof of Eriedman’s second eigenvalue Theorem and its extension to 
random lifts. Preprint, arXiv:1502.04482, 2015. 

[8] C. Bordenave and A. Guionnet. Localization and delocalization of eigenvectors for heavy¬ 
tailed random matrices. Probab. Theory Related Fields, 157(3-4) :885-953, 2013. 

[9] C. Bordenave and M. Lelarge. Resolvent of large random graphs. Random Structures 
Algorithms, 37(3):332~352, 2010. 

[10] P. Bourgade, J. Huang, and H.-T. Yau. Eigenvector statistics of sparse random matrices, 
2016. Preprint, arXiv: 1609.09022. 

[11] P. Bourgade and H.-T. Yau. The eigenvector moment flow and local quantum unique 
ergodicity. Comm. Math. Phys., pages 1-48, 2016. 

[12] S. Brooks and E. Lindenstrauss. Non-localization of eigenfunctions on large regular graphs. 
Israel J. Math., 193(1):1-14, 2013. 

[13] C. Cacciapuoti, A. Maltsev, and B. Schlein. Bounds for the Stieltjes transform and the 
density of states of Wigner matrices. Probab. Theory Related Fields, Nov 2014. 


51 


[14] N.A. Cook. On the singularity of adjacency matrices for random regular digraphs, 2014. 
Preprint, arXiv:1411.0243. 

[15] C. Cooper, M. Dyer, and C. Greenhill. Sampling regular graphs and a peer-to-peer network. 
Combin. Probab. Comput., 16(4):557-593, 2007. 

[16] I. Dumitriu, T. Johnson, S. Pal, and E. Paquette. Functional limit theorems for random 
regular graphs. Probab. Theory Related Fields, 156(3-4) :921-975, 2013. 

[17] I. Dumitriu and S. Pal. Sparse regular random graphs: spectral density and eigenvectors. 
Ann. Probab., 40(5):2197~2235, 2012. 

[18] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Spectral statistics of Erdos-Renyi Graphs 
II: Eigenvalue spacing and the extreme eigenvalues. Comm. Math. Phys., 314(3):587-640, 
2012. 

[19] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. The local semicircle law for a general class 
of random matrices. Electron. J. Probab., 18:no. 59, 58, 2013. 

[20] L. Erdos, A. Knowles, H.-T. Yau, and J. Yin. Spectral statistics of Erdos-Renyi graphs I: 
Local semicircle law. Ann. Probab., 41(3B):2279-2375, 2013. 

[21] L. Erdos, B. Schlein, and H.-T. Yau. Local semicircle law and complete delocalization for 
Wigner random matrices. Comm. Math. Phys., 287(2):641-655, 2009. 

[22] L. Erdos and H.-T. Yau. Universality of local spectral statistics of random matrices. Bull. 
Amer. Math. Soc. (N.S.), 49(3):377-414, 2012. 

[23] L. Erdos, H.-T. Yau, and J. Yin. Bulk universality for generalized Wigner matrices. Probab. 
Theory Related Fields, 154(l-2):341-407, 2012. 

[24] J. Friedman. A proof of Alon’s second eigenvalue conjecture and related problems. Mem. 
Amer. Math. Soc., 195(910):viii-|-100, 2008. 

[25] J. Friedman, J. Kahn, and E. Szemeredi. On the second eigenvalue of random regular 
graphs. In Proceedings of the Twenty-first Annual ACM Symposium on Theory of Comput¬ 
ing, STOC ’89, pages 587-598, New York, NY, USA, 1989. ACM. 

[26] L. Geisinger. Convergence of the density of states and delocalization of eigenvectors on 
random regular graphs, 2013. Preprint, arXiv: 1305.1039. 

[27] J. Huang, B. Landon, and H.-T. Yau. Bulk universality of sparse random matrices, 2015. 
Preprint, arXiv: 1504.05170. 

[28] D. Jakobson, S.D. Miller, 1. Rivin, and Z. Rudnick. Eigenvalue spacings for regular graphs. 
In Emerging applications of number theory (Minneapolis, MN, 1996), volume 109 of IMA 
Vol. Math. AppL, pages 317-327. Springer, New York, 1999. 

[29] W.B. Johnson, G. Schechtman, and J. Zinn. Best constants in moment inequalities for 
linear combinations of independent and exchangeable random variables. Ann. Probab., 
13(l):234-253, 1985. 


52 


[30] H. Kesten. Symmetric random walks on groups. Trans. Amer. Math. Soc., 92:336-354, 
1959. 

[31] A. Knowles and J. Yin. Eigenvector distribution of Wigner matrices. Probab. Theory 
Related Fields, 155(3-4):543-582, 2013. 

[32] A. Knowles and J. Yin. The isotropic semicircle law and deformation of Wigner matrices. 
Comm. Pure Appl. Math., 66(11):1663-1750, 2013. 

[33] A. Knowles and J. Yin. Anisotropic local laws for random matrices, 2014. Preprint, 
arXiv:1410.3516. 

[34] B. Landon and H.-T. Yau. Convergence of local statistics of Dyson Brownian motion, 2015. 
Preprint, arXiv: 1504.03605. 

[35] B.D. McKay. The expected eigenvalue distribution of a large regular graph. Linear Algebra 
Appl, 40:203-216, 1981. 

[36] B.D. McKay. Asymptotics for symmetric 0-1 matrices with prescribed row sums. Ars 
Combin., 19(A):15-25, 1985. 

[37] S.J. Miller and T. Novikoff. The distribution of the largest nontrivial eigenvalues in families 
of random regular graphs. Experiment. Math., 17(2):231-244, 2008. 

[38] I. Oren and U. Smilansky. Trace formulas and spectral statistics for discrete Laplacians on 
regular graphs (II). J. Phys. A, 43(22):225205, 13, 2010. 

[39] D. Puder. Expansion of random graphs: new proofs, new results. Invent. Math., 201(3):845- 
908, 2015. 

[40] Z. Rudnick and P. Sarnak. The behaviour of eigenstates of arithmetic hyperbolic manifolds. 
Comm. Math. Phys., 161(1): 195-213, 1994. 

[41] P. Sarnak. Recent progress on the quantum unique ergodicity conjecture. Bull. Amer. 
Math. Soc. (N.S.), 48(2):211-228, 2011. 

[42] U. Smilansky. Discrete graphs—a paradigm model for quantum chaos. In Chaos, volume 66 
of Prog. Math. Phys., pages 97-124. Birkhauser/Springer, Basel, 2013. 

[43] T. Tao and V. Vu. Random matrices: universal properties of eigenvectors. Random Matrices 
Theory Appl, 1(1):1150001, 27, 2012. 

[44] L.V. Tran, V.H. Vu, and K. Wang. Sparse random graphs: eigenvalues and eigenvectors. 
Random Structures Algorithms, 42(1):110-134, 2013. 

[45] N.C. Wormald. Models of random regular graphs. In Surveys in combinatorics, 1999 (Can¬ 
terbury), volume 267 of London Math. Soc. Lecture Note Ser., pages 239-298. Cambridge 
Univ. Press, Cambridge, 1999. 

[46] S. Zelditch. Uniform distribution of eigenfunctions on compact hyperbolic surfaces. Duke 
Math. J., 55(4):919-941, 1987. 


53 


